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We determine baryon-baryon bound states in 3 + 1 dimensional SU(3) lattice QCD with two 
flavors, 4x4 spin matrices, and in an imaginary time formulation. For small hopping parameter, 
k > 0, and large glueball mass (strong coupling), we show the existence of three-quark isospin 1/2 
particles (proton and neutron) and isospin 3/2 baryons (delta particles), with asymptotic masses 
—3 In ft and isolated dispersion curves. We only consider the existence of bound states of total isospin 
/ = 0, 3. Using a ladder approximation to a lattice Bethe-Salpeter equation, baryon-baryon bound 
states are found in these two sectors, with asymptotic masses —6 In k and binding energies of order 
k 2 . The dominant baryon-baryon interaction is an energy-independent spatial range-one potential 
with an 0(k, 2 ) strength. There is also attraction arising from gauge field correlations associated with 
six overlapping bonds, but it is counterbalanced by Pauli repulsion to give a vanishing zero-range 
potential. The overall range-one potential results from a quark, antiquark exchange with no meson 
exchange interpretation; the repulsive or attractive nature of the interaction does depend on the 
isospin and spin of the two-baryon state. 
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1. INTRODUCTION AND RESULTS 

Our aim in this paper is to obtain the low-lying energy- 
momentum (EM) spectrum such as particles and their 
bound states in SU(3) lattice quantum chromodynamics 

(QCD) [JHHEl. We treat a lattice Q CD model in 
3+1 dimensions, with 4x4 Dirac spin matrices, and 

with two-flavor quarks (up and down). This model has 
a global SU(2) flavor (isospin) symmetry, and accommo- 
dates twenty, three-quark, one-particle states (baryons), 
and their antiparticles, associated with the proton (p), 
the neutron (n) and the delta (A) particles. 

This work is an important part of a recent series of 
papers given in Refs. [J IE IS 0, IE El E] where, treat- 
ing less complex QCD models, we started a research 
program aiming at developing a building up process 
and techniques to understand, from first principles, the 
hadronic particles and their bound states, bridging the 
gap between QCD and nuclear physics. The existence of 
baryons and mesons is a manifestation of confinement. 
Our main goal is to understand the mechanisms respon- 
sible for attraction and repulsion between the hadronic 
particles, when and how bound states occur and how 
their binding is related to the effective Yukawa meson- 
exchange theory. 

Here, as in previous papers (see Refs. imam 

IT(il|). the analysis is performed in an imaginary-time for- 
mulation of lattice QCD, and is restricted to the strong 
coupling regime (small hopping parameter < k <C 1, 
and large glueball mass f3 = l/gfi <C ft). Even though 
the strong coupling regime is far from the scaling limit 
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(continuum quantum field theory), these models do ex- 
hibit particles and we hope, as for confinement, that the 
main features of attraction and repulsion of the contin- 
uum model are present in the strong coupling regime. 

Of course, there is the pointwise repulsive interaction 
originating from Pauli exclusion. As expected, there is 
also a nonlocal interaction arising from a quark-antiquark 
(not a meson particle!) exchange mechanism like in the 
effective Yukawa theory. However, the attractive or re- 
pulsive nature of this interaction does depend on the spin 
and isospin of the two-baryon state. How this depen- 
dence occurs is to be understood. Not less important, our 
analysis displays another essential mechanism of attrac- 
tion arising from gauge field correlation effects associated 
with integration over (six) overlapping gauge bonds: if 
these gauge correlations are not present the bound states 
cease to exist. 

Our general approach to detect particles and bound 
states in the low-lying part of the spectrum is to de- 
rive spectral representations for the two- and four-baryon 
functions which allow us to relate their complex momen- 
tum singularities to the EM spectrum via a Feynman-Kac 
(F-K) formula. For other approaches to the determina- 
tion of the spectrum both at the theoretical and numer- 
ical levels, see RefsJl|Jn|Jl|JJl E E E E E E 
EE1 and Refs. @ IH H dealing with numerical 
simulations. 

It must be pointed out that the detection of particle 
(and their bound states) masses through the exponential 
decay rate of suitable correlations, without a spectral rep- 
resentation, is meaningless, and the resulting values may 
be far from the correct ones, especially in cases where 
degeneracies are broken with small spectral separations. 

To determine the bound state spectrum, in a first step, 
we use a ladder approximation to a lattice version of the 
Bethe-Salpeter (B-S) equation for suitable four-baryon 
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correlations. However, our work does incorporate the 
main ingredients to go beyond this ladder approximation 
and analytically validate our results to the full model 
in the strong coupling regime, similar to what is done 
in some classical and stochastic spin systems (see Refs. 

heed 

In our previous papers, the simplest case in which a 
two-baryon bound state appeared is in the two-flavor to- 
tal isospin / = 0, 1 sectors with 2x2 spin matrices in 
2 + 1 dimensions (see Ref. 0); there are no bound states 
for 7 = 2,3. However, this model is not complex enough 
to accommodate states like protons and neutrons in the 
one-particle spectrum. 

Here, considering a more realistic lattice QCD model, 
it is shown that there are twenty one-baryon states with 
asymptotic masses given by —3 In n — 3k 3 /4 + 0(k 6 ) and 
isolated dispersion curves (the upper gap property). The 
upper gap property is unknown in the Hamiltonian for- 
mulation of the model. 

Concerning the two-particle spectrum, we determine 
the two-baryon bound states for the total isospin 1 = 
and 1 = 3 sectors and below the two-baryon threshold, 
which is given by twice the smallest of the baryon mass. 
We show that there are diverse two-baryon bound states. 
Their asymptotic masses are — 61nK and their binding 
energies are of order k 2 . In the 1 = sector, the most 
strongly bound, two-baryon bound states correspond to 
a superposition of p — n and A — A total spin 5=1 
states, and also A — A, 5 = 3 states. The more weakly 
bound, bound states are associated with a superposition 
of p — n and A — A total spin S = states and, in 
addition, A — A, S = 2 states. In contrast to the 1 = 
states, we find that for the maximum isospin 1 = 3 sector 
there are also strongly bound and weakly bound bound 
states in the lowest spin sectors S = 0, 1 and no bound 
states if S = 2, 3. These results are in agreement with our 
previous results of Ref. 9] that the attraction between 
the two particles decreases with increasing /. 

It is important to note that although our results are 
obtained using fairly complicated machinery, in the end, 
a simple picture emerges for the formation of baryon- 
baryon bound states. There is a correspondence between 
the lattice B-S equation, in relative coordinates, and a 
lattice resolvent operator Schrodinger equation. The two- 
baryon system in relative coordinates behaves approxi- 
mately like that of a non-relativistic one-particle lattice 
Hamiltonian T + V with lattice kinetic energy — k 3 A^, 
where A^ is the spatial lattice Laplacian and the poten- 
tial energy V is k 2 V' , the quasi-meson exchange space 
range-one potential which dominates the kinetic energy 
for small k. The attractive or repulsive nature of the in- 
teraction depends on the isospin, spin spectral structure 
of V at a single site of space-range one. Because of the 
k 3 dependence of the kinetic energy T and the k 2 depen- 
dence of the potential energy, there is no minimal critical 
value of the interaction strength needed for the presence 
of a bound state. This is in contrast to the case of the 
Birman-Schwinger (see Ref. j2^]) bound for a — Ag + W. 



If V is everywhere positive, there is a critical negative A 
needed for a bound state to occur in three space dimen- 
sions. 

A brief report on the results given in this paper ap- 
pears in Ref. [29|. The present article is meant to be 
much more self-contained, describing in a relatively com- 
plete way the details about our methods and techniques. 
It is organized as follows. In Section [21 we mathemati- 
cally define the model, describe the standard construc- 
tion of the underlying physical Hilbert state 7i, and in- 
troduce the fields which naturally emerge as those quan- 
tities associated with the physical baryon particles. To 
do this, we use a hyperplane decoupling expansion as 
done for quantum field theories in the continuum (see 
[3fj|). but adapted to the lattice. In Section |3 the one- 
particle analysis is done, and the baryon-baryon bound 
states are determined below the two-baryon threshold in 
Section 21 In order to keep the text as fluid as possible, 
trying to separate the main ideas of our methods from 
unavoidable sometimes heavy technical details, Sections 
13 and 0] are subdivided according to their principal log- 
ical steps and details about specific points are given in 
the Appendices A-F. The similarity between the lattice 
B-S, in relative coordinates, and the lattice resolvent op- 
erator Schrodinger equation is exploited in Section^ Fi- 
nally, in Scction|51 we end the paper with conclusions and 
prospective ideas, and a list of important open problems. 

To close the Introduction, it is useful to make some 
comments on the seemingly awkward form of the two- 
and four- baryon correlation functions that we use in the 
sequel and which play an important role in our method. 
We explain the guiding principle for obtaining the form of 
correlation functions, which also holds in the continuum 
spacetime. Consider the physical Hilbert space matrix 
elements (Xi e where H is the energy operator, 

e -n\t\ j g imaginary-time evolution operator (semi- 
group) and x an d V* are both in the subspace of one- 
or two-baryon states. We have a F-K formula at our 
disposal and applying it to the two time orderings in t 
gives the form of the correlation functions. It is these 
correlation functions, that have appropriate imaginary- 
time decay properties, which allows our analysis to go 
through. 



2. THE MODEL, THE PHYSICAL HILBERT 
SPACE AND BARYON FIELDS 

We now introduce our SU(3), d + 1 dimensional, 
imaginary-time lattice QCD model, where d = 3 is the 
space dimension. The ^partition function is given for- 
mally by Z = J e ~ sw ''^ 9 ) dipdtpd(x(g) , and for a func- 
tion F(il>,ip,g), the normalized correlations are denoted 
by 

(F$,1>,g)) = ~ / F{^,g)e- s ^U^d^dii{g). 

(2.1) 
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The model action S = S(tp, ip, g) is Wilson's action 

S = - ^ $a, a ,f(u)T° e g {g U M+ae»)ab^b,{3j{u + ffe M ) 



x(g P ), 



(2.2) 

where, besides the sum over repeated indices a, f3 = 
1,2,3,4 (spin), a = 1,2,3 (color) and / = +1/2,-1/2 = 
+,— (isospin), the first sum runs over u — (it , it) = 
(wVW 3 ) eZ*E {±1/2, ±3/2, ±5/2...} x Z 3 , a = 
±1 and \i = 0,1,2,3. Here, we are adopting the la- 
bel for the time direction and e M , fj, = 0,1,2,3, de- 
notes the unit lattice vector for the /i-direction. The 
choice of the shifted lattice for the time direction, avoid- 
ing the zero-time coordinate, is so that, in the contin- 
uum limit, two-sided equal time limits of quark Fermi 
fields correlations can be accommodated. At each site 
u E Zq, there are fermionic Grassmann fields ipaaf(u), 
associated with a quark, and fields ip aa f(u), associated 
with an anti-quark, which carry a Dirac spin index a, an 
SU(3) color index a and isospin /. We refer to a = 1, 2 
as upper spin indices and a = 3, 4 or + or — respec- 
tively, as lower ones. For each nearest neighbor ori- 
ented lattice bond < u,u± e' 1 > there is an SU(3) ma- 
trix U {g u ,u±e^ ) parametrized by the gauge group element 
g u ,u±e» and satisfying Uigu^+e^y 1 = U(g u+e » >u ). Asso- 
ciated with each lattice oriented plaquette p there is a pla- 
quette variable x(U(g p )) where U(g p ) is the orientation- 
ordered product of matrices of SU(3) of the plaquette 
oriented bonds, and x is t ne trace. For notational sim- 
plicity, we sometimes drop U from U(g). Concerning 
the parameters, we take the quark-gauge coupling or 
hopping parameter k > 0. Also, go > describes the 
pure gauge strength and M = M(m, n) = (m + 2n)I^ 1 
I n being the n x n identity matrix. Given k, for sim- 
plicity and without loss of generality, m > is cho- 
sen such that M a p = 8 a p, meaning that m + 2k = 1. 
Also, we take T ±e ^ = —I4, ± 7 M , where 7^, satisfy- 
ing {7^,7,,} = 26^1^, are the 4x4 hermitian trace- 

h 
^ -h 

ia ' ' 
-ia j 

the hermitian traceless anti-commuting 2x2 Pauli spin 
matrices. The measure dfi(g) is the product measure 
over non-oriented bonds of normalized SU(3) Haar mea- 
sures (see 31]). There is only one integration variable 
per bond, so that g uv and g~^ are not treated as dis- 
tinct integration variables. The inte gra ls over Grass- 
mann fields are defined according to |32| . For a poly- 
nomial in the Grassmann variables with coefficients de- 
pending on the gauge variables, the fermionic integral is 
defined as the coefficient of the monomial of maximum 
degree, i.e. of \[ ul ipe(u)ip e (u), I = (a, a,/). In Eq. 
(|2.1|) . dtp dip means Y[ u 1 dipt{u) dipi(u) such that, with 
a normalization A/"i = (1), we have (ip£ 1 (x) tpe 2 (y)) = 



less anti-commuting Dirac matrices 7 = 



and j J 



and a\ j = 1,2,3, denotes 



(1/A/i) / ipM^tM^^'^ *' 3( " )0 ' 3,< ^ w ## 
= a2 6 aia2 S f 1 f 2 S(x — y), with a Kronecker delta for 
space-time coordinates, and where O is diagonal in the 
color and isospin indices. 

Throughout the paper, we work in the strong coupling 
regime i.e. we take < t/ 1 < k < 1, m = 1 - 2k < 1. 
With this restriction on the parameters, and since we 
have chosen r ±e matrices within the two-parameter 
family described in Ref. [13j, positivity is preserved and 
there is a quantum mechanical Hilbert space H. of physi- 
cal states (see below) . Furthermore, the condition m > 
guarantees that the one-particle free Fermion dispersion 
curve increases in each positive momentum component 
and is convex for small momenta. 

The action of Eq. 12.21 is invariant by the gauge trans- 
formations given by, for x € Z^ and h(x) € SU(3), 



U(g 



x-\-e^ ,x . 



ip{x) 1— > h(x)ip(x) 
ip(x) 1— > ip(x) [ft.(x)] _1 
i-> [h(x + e li )]~ 1 U(g x+e ^ iX )h(x) 



(2.3) 



Concerning the global SU(2) isospin symmetry of the ac- 
tion, we follow the treatment given in Ref. ,10]. We 
treat isospin symmetry at the level of correlation func- 
tions (rather than as operators in the physical Hilbert 
space) in Appendix A. Other symmetries of the action 
in Eq. (|2.2|l . such as time reversal, charge conjugation, 
parity, coordinate reflections and spatial rotations, which 
can be implemented by unitary (anti-unitary for time re- 
versal) operators on the physical Hilbert space H, were 
treated in Ref. |5| and do not affect the isospin indices. 
For k — 0, we recover the spin SU(2) symmetry, sepa- 
rately in the lower and the upper spin indices. Further- 
more, in Appendix B, a new time reflection symmetry, 
to be distinguished from time reversal, of the correlation 
functions is shown to be useful for obtaining additional 
relations between correlations. 

By polymer expansion methods (see Refs. j0,H3|), the 
thermodynamic limit of correlations exists and truncated 
correlations have exponential tree decay. The limiting 
correlation functions are lattice translational invariant. 
Furthermore, the correlation functions extend to analytic 
functions in the coupling parameters k and (3 = 1 /<7n- For 
the formal hopping parameter expansion, see Refs. [2J,|3|. 

The underlying quantum mechanical physical Hilbert 
space Tt and the EM operators H and P J , j = 1, 2, 3 are 
defined as in Refs. We start from gauge invariant 

correlations, with support restricted to u° — 1/2 and 
we let Tf°, Tf , i = 1,2,3, denote translation of the 
functions of Grassmann and gauge variables by a; > 0, 
x — (2; 1 , x 2 , a; 3 ) £ 1? . For F and G only depending on 
coordinates with u° = 1/2, we have the F-K formula 



(G,TfffffffF) n 



[Tf T s F]OG) , (2.4) 



where T x = Tf Tf Tf and O is an anti-linear operator 
which involves time reflection. Following Ref. [13, the 
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action of on single fields is given by 

<dlpaa(u) = (jo)ai3lpal3(tu) , 

where t(u ,u) = (— u°,u), for A and B monomials, 
Q(AB) = Q(B)Q(A); and for a function of the gauge 
fields Of ({g uv }) = f*({g(tu)(tv)}), u,v e Zf+\ where 
* means complex conjugate. 9 extends anti-linearly to 
the algebra. For simplicity, we do not distinguish be- 
tween Grassmann, gauge variables and their associated 
Hilbert space vectors in our notation. As linear opera- 
tors in Ti, T M , n = 0, 1, 2, 3, are mutually commuting; To 
is self-adjoint, with — 1 < To < 1, and Tj = \.2^ are uni- 
tary. So, we write fj = e lP ' and P = (P 1 , P 2 , P 3 ) is the 
self-adjoint momentum operator. Its spectral points are 
p £ T 3 = (— 7r, 7r] 3 . Since Tg > 0, the energy operator 
H > can be defined by Tq = e~ 2H . We call a point 
in the EM spectrum associated with spatial momentum 
p = as a mass and, to be used below, we let £ (A , A) be 
the product of the spectral families for the operators Tn, 
Ti, T2 and T3. By the spectral theorem (see Ref. [28j . 
we have 



T 



X°dE (X ) 



T 



i=l,2,3 



X°dFj(X j ) 



so that £(A°, A) = E (X a ) J]? F 3 (X^>). The positivity con- 
dition (FQF) > is established in but there may 
be nonzero F's such that (FQF) = 0. If the collection 
of such F's is denoted by AT, a pre-Hilbert space Ti' can 
be constructed from the inner product (GQF) and the 
physical Hilbert space Ti is the completion of the quo- 
tient space Ti'/fif, including also the cartesian product 
of the inner space sectors, the color space C 3 , the spin 
space C 4 and the isospin space C 2 . 

Here, we analyze the one and two-baryon sectors of Ti. 
Points in the EM spectrum are detected as singularities 
in momentum space spectral representations of suitable 
two- and four-baryon field correlations given below. 



3. ONE-PARTICLE SPECTRUM 

In this section, we restrict our attention to the sub- 
space Ti C Ti generated by an odd number of ^'s where 
the baryons lie (the tilde means the presence or absence 
of a bar). The mesons lie in the even subspace Ti e - We 
treat the one-baryon states adapting the methods of Refs. 
5,9]. 



3.1. Baryon fields and particle identifications 

The decoupling of hyperplane method (see Refs. 0,0) 
reveals that the gauge invariant barred fields listed be- 
low create baryon particles. We point out that the same 



method applies to QCD models with any number of fla- 
vors. In particular, for the SU(3) invariant three flavor 
case, the method would give the eightfold way. Here, the 
baryon particles are labeled by one-baryon total isospin 
(I s * index), the one baryon z-component isospin (2 nd in- 
dex) and z-component of total spin (3 rd index). 

Let a hat " mean the presence or absence of a bar si- 
multaneously for all ^'s. The decoupling of hyperplane 
expansion shows that the baryon fields are given by (with 
all fields evaluated at the same point) 



2 2 ~T 

Bi -1 ±1 



B311 = 

2 2 2 

B 3 i-i 

2 2 2 

B3 1 ±3 
2 2 2 



(3.5) 



= §^§ ('<Pa± + 1pb T - ~ i>a±-lpb T +)ll>c±+, 

^(ip a++ 4>b — +2?/> a+ _?/> fc _ + )v3 c++ , 

= ^(2$ 0++ &__-^ 0+ _&_+)$ c -+, (3.6) 



B 3 

2 


-11 = 
2 2 


= ^t(^a++^b— +V-V+- 


-i>b— 




B, 

2 


-1 —1 

2 2 




—i>b- 


-+)lpc— 


B, 

2 


-1 ±3 
2 2 


= ^7=^a±--0b±-^ C ±+, 






B 3 

2 


-3 ±1 
2 2 


= ^=4>a±-Tpb±--lpc T - 






B 3 

2 


-3 ±3 
2 2 


= ^ tp a ±-^b±~A±- 






B, 

2 


3 ±1 = 
2 2 


= §^§ ^a± + A± + Tjjc T + , 






B 3 

2 


3 ±3 = 
2 2 


= ^a± + 1pb± + 1pc±+ , 







(3.8) 
(3.9) 



where the index a in tp at a,f is a lower spin component 
a = 3,4= +,-. 

The baryon fields in Eqs. H3.5fH3.9[l obey the normal- 
ization (BiBj)^ = —Sij, where i,j are collective indices 
and the superscript ^ denotes k = in the hopping term 
in the action S given in Eq. 11'. 21) . Below, the superscript 
means the coefficient of c", n = 0, 1, ... We recall that, 
at k = 0, our model also has a global SU(2) spin invari- 
ance in the space of lower and upper components. 

A brief description of the hyperplane decoupling ex- 
pansion method, and an important product structure 
for special third order k derivatives, at n = 0, for the 
two-baryon correlation function defined below is given in 
Appendix B [see Eq. (|B.4|l ]. together with the isospin 
identification for the above baryon fields with physical 
particles. The method reveals the fields and associated 
particles so that no guesswork or ansatz has to be used. 

In analogy with the two-flavor SU(3) gauge QCD in 
the continuum with only the up (electric charge +2/3, z- 
component isospin +1/2) and the down (electric charge 
— 1/3, z-component isospin —1/2) quarks, we identify the 
particles associated with the barred fields in Eqs. (|3.5|) - 
(|3.9() . The proton (p±) and the neutron (n±) are, re- 
spectively, described in Eq. 1)3.5(1 . A + , A and A - are 
identified with the fields in Eqs. H3.6H3.8fl . Finally, 
A ++ corresponds to the baryon fields in Eq. 1)3.9(1 . The 
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corresponding z-component isospin conjugated particles 
are described by the unbarred fields. 

We remark that, although not obvious, the fields in 
Eqs. I|3.5[l - i|3.9[) are the same as those obtained by the 
construction of gauge-invariant (colorless) wave functions 
as products of spin symmetric, flavor symmetric wave 
functions, except for the proton and neutron. For the 
proton and neutron, wave functions of Eq. (|3.5|) are the 
same as the sum of products of wave functions that are 
antisymmetric in both spin and flavor in the pairs 12, 13 
and 23. 



3.2. Two-baryon correlations and the one-particle 
analysis 

The baryon-baryon correlation function is defined by 
(x here denotes the characteristic function) 

G M2 (u,v)= (B tl (u)B l2 (v)) Xu"<v" 

- (B ei (u)B e2 («))* Xu°>v° = G hi2 (u-v), 

_ (3.10) 

where £ = (I, I z , s) is a collective index. By isospin sym- 
metry Gi x i 2 (x — u — v) is diagonal in /, I z and, for / 
fixed, independent of I z . The consequences of isospin 
symmetry for the two-point functions, as well as for the 
four-point functions given below, are discussed in Ap- 
pendix A, where important orthogonality relations and 
correlation function identities are derived. Also, the z- 
component of spin s is restricted to ±1/2, for 7=1/2 (for 
p and n) or ±3/2, ±1/2, for I = 3/2 (the A particles). 

Setting x° — u°—v° ^ 0, we have the following spectral 
representation for G, with Bi = B^(l/2, 0), 



Ge 1 e 2 (x) = - 



r^o\\x°\-i e iX.x 
x d\(Be 1 ,£(X°, \)Bi 2 )h ; 



(3.11) 



for x € Z d+1 , x° 0, and is an even function of x. 

Following Proposition 1 of Ref . [f| , we verify how Eq. 
(|3.11() is obtained. First, we derive a spectral represen- 
tation for u° < v°. With this time ordering, we have 

Ge x e 2 (u, v) = (B £l [u)B l2 («)> 

= -([T " - u °- 1 f ff - s ^ 2 (l/2 ) 0)}B ei (-1/2,0)). 

But B £ (-l/2,0) = 9^(1/2,0), so that 

G tlt Jiu,v)^-([TS - u0 - 1 r , - a BUl/2, 0)}OB ei (1/2,0)). 

By parity symmetry, Bi(w°,0) has parity —1 so that 
v — u can be replaced by u — v here and be- 
low. The F-K formula of Eq. (J23J) gives the 
result. For vP > v , we write G^iJujV) = 

-([T tl °-'' O - 1 f a -^ 1 (-l/2,0)]eB, 2 (l/2,0))* and use 
the F-K formula. For u° = v° there is no spectral rep- 
resentation. By going backwards, the correlation of Eq. 
(|3.10|) is obtained from the Hilbert space inner product. 



With this, upon separating the time zero contribution, 
the Fourier transform G^i 2 (p) = X^ez 4 G^i 2 (x)e~ lp - x 
of G^^ix) admits the spectral representation 

Ge lta (p) = G eie M-(2n) 3 f\f(p ,\ )d x «a fJie2 (\ ), 

(3.12) 

where 



T 3 



S(p-X)d x od x (B ei ,£(X°,X)B i2 ) 



H 



(3.13) 

y) , and we set 



with f(x,y) = (e lx - y) 1 + (e 
G(P) = Eg e-^- s G(x° =0,x). 

The spectral representation given in Eq. I|3.12|l allows 
us to identify complex p singularities of G(p) with points 
in the EM spectrum. 

Using symmetries, Gg^^p) is real diagonal and inde- 
pendent of isospin and color indices, diagonal and inde- 
pendent of the spin, at least up to order k 6 . 

From isospin symmetry, we see that the 20 x 20 two- 
point function matrix decomposes into blocks. There is 
one 2x2 block for the proton p and the same for the 
neutron n, and there are four 4x4 blocks associated to 
the A's. Turning to the determination of the dispersion 
curves, as G^ 2 is diagonal in / and I z , and for fixed /, 
independent of I z , we consider / and I z fixed and denote 
the spin dependent two-point function by G SlS2 . The dis- 
persion curves w(p) = w(p, k) are given by the zeroes of 
the determinant of T SlS2 (p), for p° in the positive imagi- 
nary axis, i.e. 

detf SlS2 (p° = iw(p),p) = 0. 

In Appendix B, we show that the 20 x 20 two-point 
function matrix T SlS2 (p) is not diagonal at p ^ (it is 
diagonal for p = 0!). Thus, if we want to obtain the dis- 
persion curves, this property prevents us to follow closely 
the analysis given in Ref. [5| for the mass spectrum. Us- 
ing a different method, we show that there are 2 (4) dis- 
persion curves, not necessarily distinct, for 1 = 1/2 (3/2), 
denoted by Wj(p), j = 1, 2 for / = 1/2 and j = 1, 2, 3, 4 
for / = 3/2; the masses rrij are given by rrij = Wj(p = 0). 
The twenty approximate one-particle dispersion curves 
are identical up to and including 0(k 5 ) and, with 



»2 = 2^(l-cosp l ), 



(3.14) 



are given by 



Wj(p) = Wj(p, n) = [-3hiK - 3k 3 /4 +p 2 e K, 3 /8] + r 3 (p) , 

(3.15) 

where rj(p) is of C(k 6 ), and its rj(p = 0) at 0(k 6 is 
zero. Thus, if there is a mass splitting, it is at least 
of 0(k 7 ). As also shown in Appendix B, for the case 
of 2 x 2 blocks associated with the proton and neutron, 
since the 2x2 matrices are diagonal by a parity-charge 
conjugation-time reflection symmetry, we can apply the 
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more powerful method of Ref. p|, which includes an 
application of the analytic implicit function theorem, to 
show that the corresponding (identical) rj (p) are analytic 
in k and in each component p l=1 > 2 * 3 of the momentum p. 

The spectral measure dxoap )SlS2 (A ) of Eq. (|3.13|) has 
the decomposition 

dxoa p , slS2 (\°) = £\ [Zj tSlS2 (p)S(X° - e- w ^)]d\° 
+d\ov SlS2 {\°,p) , 

where the first term corresponds to the one-particle con- 
tributions, and 

% SlS2 (p)- 1 = -(27r) 3 e^®-^(/ = ix,P)\ x = Wi (f) ■ 

Let r siS2 (it, v) denote the convolution inverse of the two- 
point function G SlS2 (u,v). We can show that T SlS2 (u, v) 
decays faster than G SlS2 (u,v), and hence the Fourier 
transform f SlS2 (p) = [G]j^ S2 (p) of T SlS2 (x = u - v) = 
r siS2 (it, v) has a larger region of analyticity in p°. Thus, 
the cofactor [cof r] SlS2 (p)/det[r(p)] = Tj i 1 S2 (p) provides 

a meromorphic extension of G SlS2 (p) up to near the two- 
particle threshold (precisely, the meson-baryon thresh- 
old, which is sa — 51n«;), and the one-particle EM spec- 
trum occurs, for each p, as points given by the p° imag- 
inary axis zeroes of det[F(p)]. The two-point function 
convolution inverse T(x) is defined by 



and then 



(1 + Gj 1 G n ) 1 Gj 1 . 



(3.16) 



using a Neumann series, where Gd is the diagonal part 
of G given by 



Gd,s lS2 (u,v) = G SlS2 (u,u)5 SlS2 S u 



(3.17) 



and G n is the remainder 



G n ,s lS2 (u, v) = G SlS2 (u, v) - G d , Sl s 2 (u, v) . (3.18) 

By the decoupling of hyperplane method, we show the 
decay of G and T; G^ 1 , G n and T are bounded as ma- 
trix operators on ^(C 20 x Z^); for small k, G n is small 
in norm so that the Neumann series in Eq. (|3.16|) con- 
verges. To obtain the results on the dispersion curves 
Wj(p), up to order k 6 , we need detailed information on 
the small distance behavior of r siS2 , which are obtained 
by showing explicit cancelations in the Neumann series 
and go beyond the hyperplane decoupling method (see 
Ref. [|| and Appendix B for more details). 

We now give an intuitive argument for obtaining the 
mass and dispersion curves. The precise value requires 
a more refined analysis adapting the method used in 
Ref. p|, and depends on the short distance behavior of 
r siS2 (x), which in turn depends on that of G SlS2 (x), via 
the Neumann series. By explicit calculation, G SlS2 (0) = 

" 3, j = 1,2,3, 
) such that 



-1, G SlS2 (±e°) = -k 3 , G SlS2 (±eJ 
(correcting the misprinted values in Ref. 



5] 



G SlS2 (p) = [-l-2cosA 3 -^r^ cosp>]6 SlS2 +0( K 6 ) , 



f SlS2 (p) = [-1 + 2cos/k 3 + — cosp>}5 SlS2 +0(n 6 ) 

3=1 



Setting r siS2 (p) = for p° — iwj(p), we obtain the dis- 
persion curves of Eq. I|3.15|l . which are shown to be the 
correct ones up to C(k 5 ) by adapting the methods of Ref. 

a 

Finally, with this, the A support of the measure 
dv SlS2 (\°, A) is contained in |A°| < |k| 5_c . 
The decomposition of G SlS2 (p) is given by 



G SlS2 (p) = G SlS2 (p) - J2j I z j,s iS2 (p) 
1 



1 



e -ip° _ e wj(p) 

x d x ov^ SlS2 (X°) 



3 f 1 f(„Q \0\ 



(2 7 r) 3 /_ 1 /(p°,A°) 

(3.19) 



4. TWO-BARYON SPECTRUM 

In this section, we restrict our attention to the sub- 
space He generated by an even number of ^'s. We treat 
the two-baryon states adapting the methods of Ref. 
(see also 0). 



4.1. Two-baryon states and four-baryon 
correlations in the individual and total spin basis 

To determine baryon-baryon bound states, we con- 
sider the subspace of the subspace TL e c H generated 
by the product of the twenty one-baryon fields given 
by Bi 1 (xi)Bi 2 (x2)- If we ignore possible linear depen- 
dencies, this space has dimension 400, which -upon us- 
ing isospin symmetry- decomposes into two-baryon total 
isospin / = 0,1,2,3 sectors, of dimensions 20, 108, 160 
and 112, respectively. In this subsection, to reduce the 
complexity of the analysis, here we first restrict our at- 
tention to the I — sector, where we expect to find a 
deuteron, p — n bound state. The I = 3 states will be 
treated in Section 4.4. The more algebraically complex 
7 = 1,2 sectors are presently being analyzed 0| ■ The 
isospin orthogonality relations and correlation identities 
used here are derived in Appendix A. 

The 1 = states are given by the isospin Clebsch- 
Gordan (C-G) linear combinations 



= 1/2 



t;;; 2 (x 1 ,x 2 ) = J2 c ixl B i,ii,si( x ^ B i,i 2 , s A x ^ > 

h+h=0 



3=1 



for the coupling of two isospin 1/2 baryons, and similarly, 
with (1/2 — > 3/2, Si — > tj), for the coupling of two 3/2 
isospins. We use Sj (tj) for isospin 1/2 (3/2) baryons and 
recall that the 1/2, 3/2 coupling does not give 1 = 0; and 
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we refer to this description of the two-particle states as 
the individual spin basis. 

It turns out that the dominant interaction between 
baryons admits a simpler description in terms of total 
spin. The total spin states are obtained by taking C-G 
combinations of the r's. Namely, 

Ts's 2 z {xi,x 2 ) = aiifHX{xi,x 2 ), (4.20) 

si+s 2 =Sz 

and similarly with (1/2 — > 3/2, S — > T, Sj — > U), where 
the a's are the spin C-G coefficients. Using their proper- 
ties, we find that 

Ts£(*i,**) = (-l) S+ %f z (x 2 , Xl ), 
T^ z {xi,x 2 ) = {-l) T+1 t^ z {x 2 ,x 1 ) , 

so that, for coincident points, the states vanish for 5 = 0, 
and for T — 0, 2, a kinematic statistical effect. The states 
are in the space symmetric (anti-symmetric) for 5=1, 
T = 1,3 (5 = 0, T = 0,2). 

We order the individual spin pairs as {81,82) — 
\2> 2)' v2> — 2)' ( — 2' 2)' ( — 2> — 2)' (^1^2) — (2' 2)' 

("3 1) (3 _U ,3 „3N ,1 3} ,1 _3 ) / l 3} 

V2'2'' V2' 2'' ^2' 2'' V.2'2'' ' \2' 2'' v 2'2' 1 ' 

(— |, — |), (— |, |), (— |, — |), and the total spin pairs 
as (5,5,) = (1,1), (T,T Z ) = (1,1), (5,5,) = (1,0), 
(T,T«) = (1,0), (5,5 Z ) = (1,-1), (T, T z ) = (1,-1), 
(T,T Z ) = (3,3), (3,2), (3,-3), (5,5.) = (0,0), 
{T,T Z )= (0,0), (T,T.)= (2,2), (2,1), (2,-2). From 
Eq. (|4.21() and in the total spin basis we see that the 
two baryon states are space symmetric (anti-symmetric) 
in the 1st through the 13th (last 7) components. 

The individual and the total spin basis, in the 1 = 
sector, are related by a real orthogonal transformation 
whose entries are C-G coefficients. Let A = [A] i =1 2 o 
be the 20 x 20 transformation matrix from the indi- 
vidual to the total spin basis. The nonzero elements 
of A a re: An = 1 = A45 = A75 = -4i 32 o; A27 = 
V3A0 = -4 2 i 3 = A 6 12 = A 6 1 8 ; -4 2 io = -y/2/E = 

^615! -432 = \/l/2 = -433 = -4.86 = -4sg = Ai2ie — 
-4l2 19 = Aug = -4l66 = Am = -4i9i2 = -420 16! 

-448 == v9/20 = ^417 = -A1011 = .4.1014; -44 11 = 
-yT720 = -4414; -4 97 = y/T/5 = A 913 = Ann = 
An is', Ag 10 = \/3/5 = An 15; -4io8 = \A/20 = A10 17; 

-4i43 = — yl/S = -4i69 = -4i7i3 = ^4.1918 = .42019; 

-4i5 8 = I/ 2 = -4i5i4 = -4i8 8 = -4isii and, finally, 

-4l5 11 = —1/2 = -4i5 17 = Ais 14 = -4l8 17- 

4.2. Four- baryon correlations and two-baryon 
spectrum 

As we do not know a priori which linear combination 
of states corresponds to a bound state, we consider the 
20 x 20 matrix of (truncated) four-baryon functions 

\ ( M u M 12 \ , \ 
M{xi,x 2 ;x 3 ,x A ) = I M M I {xi,X2;x 3 ,X4) , 



where, in the individual spin basis, and for 
x® = x\, x\ = x\, x\ < x®, we have 

Mn;s 1 s 2 s 3 s 4 ,(xi,X 2 ;X3,X4 : ) = (Tlis 2 {xi,X2)fl3s i (x 3 ,X4 : )), 
M 1 2- SlS2 t 3 tA x lyX 2 ;X 3 ,X i ) = (Ts{s 2 (xi,X2)f%t i {x 3 ,X 4 ,)), 
M 2 l;t 1 t 2 s 3 sJx 1 ,X2;X 3 ,X4) = {Tt[t 2 {xi,X2)flls l {x 3 ,X A )) 

and M22;t 1 t 2 t 3 t i {xi,X2;x 3 ,Xi) = { T t(t 2 {xi, x 2 ) 
X ft^t A {x 3 , Xi)); and similarly for the total spin ba- 
sis. For x\ > x 3 , Mij is defined by interchanging the 
barred and unbarred fields and taking the complex 
conjugate. As in Ref. @, Xi, X4 are now taken in Z 4 . 

has spectral representations which we use to detect 
bound states below the two-baryon threshold. We now 
obtain a spectral representation for My ;SlS2S3S4 . To sim- 
plify and make the analysis more transparent, we derive a 
spectral representation in terms of the basic baryon fields 
B^. By simply writing the suitable linear combinations 
of these representations, we obtain the spectral represen- 
tations for the four-point functions of the ti fields. 

The starting point employs the F-K formula which fur- 
nishes a 4-point correlation function in terms of Bis. For 
x° ^ 0, with x = {x° = v° - u°, x) el 4 , 

(^(1/2,^)5^(1/2,^2), (Tb)l* l-i 

x f g B i3 {l/2,u 3 )B £i {l/2,u 4 ))n = -ft lWs ! t (4 

where Gi x i 2 i 3 i A {x) = Gt l i 2 i 3 t iL {ui,U2, u 3 +x, u 4 + x) and, 
for u® = u® = u ° and W3 = u® = v , 

Ql 1 l 2 l 3 lA u l^ u 2,U 3 ,U 4 ) = 

(Bi 1 {ui)B i2 {u 2 )Be 3 {u 3 )B ti (u 4 )) Xu»<v» 

+ {Bt 1 {u 1 )Be 2 {u 2 )B e3 {u 3 )B ei {u i ))* Xu«>v° ■ 

We obtain a spectral representation for Gt^( 3 t 3 t t (x) and 
its Fourier transform Ge 1 i 2 i 3 i i {k) which allows us to re- 
late the complex k singularities of Ge. 1 £ 2 e, 3 e i {k) to the EM 
spectrum. Inserting the spectral representations for Tq 
and Ti=i,2,3, and taking the Fourier transform gives 

Gi 1 i 2 w i {k)=Gi lW Jk) - (27r) 3 Xlj T 3 f{k°,X°) 

xS{k - X)d x «d x {B ei (1/2, u^Be, (1/2, u 2 ), 
xS{X ,X)Be 3 {l/2,u 3 )B u {l/2,u 4 ))n, 

where Ge 1 e 2 e 3 e 4 { k )^J2xGT^ e ~ lk ' x Sii£ 2 e 3 £ i { x0 =0,f). The 
singularities of Ge 1 e 2 e 3 e i {k), for k = {k = ix,k = 0) and 
e ±x < 1, are points in the mass spectrum. 

To analyze Gi 1 e. 2 e. 3 e. i ik), we follow the method of anal- 
ysis for spin models as in Ref. [2(|. To use a notation 
closer to that of Ref. [2(|, we relabel the time direc- 
tion coordinates in ft 1 £ 2 £ 3 £ 4 (a;) by integer labels, with 
u { \ — 1/2 = x\, ui = Xi, i = 1,...,4, and denote by 
De 1 e 2 e 3 e 4 the shifted correlation functions. Next, we write 
De 1 e 2 e 3 e i (xi,X2,x 3 + x,x 4 + x), x\ = x\ and x% = x\, 
x° = x® — x\, where Xi and x are points on the Z 4 lat- 
tice. Now, by translation invariance, we pass to differ- 
ence coordinates and then to lattice relative coordinates 
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£ = x 2 — x±, "f] — x^ — x 3 an( i r — x 3 ~ x 2 (see Figure 
1 below). These lattice relative coordinates are a substi- 
tute for the usual center of mass and relative coordinates 
customarily used in the continuum. After doing this, we 
obtain 

Dl 1 £ 2 £ 3 e 4 (x 1 ,X 2 ,X 3 + X, Xi + x) 

= Di 1 e 2 £ 3 £ i (0,x 2 - xi,x 3 -xi + x,X4-xi + x) 
= D lll2l3ii (£,f],T + x) , 

and 

GtititaU (k) = e tfc - r De li2i3 e 4 (£, ff, k) , 

where D^e^^rf, k) = Erez* D tit 3 t 3 t*(€i V, r)e" lk - T . 
Explicitly, we have 

D e 1 e 2 e 3 e 4 (x 1 ,X2,xs, x±) = (B ei (x\ + 1/2, fi) 

x B l2 {x\ + l/2,x 2 )B i3 (x° 3 + 1/2,53) 

x Bt 4 {x\ + 1/2, Xi)) Xx°<x° 

+ (B tl (x° + 1/2, x x )B l2 {x% + 1/2, x 2 ) 

x B e3 (x° 3 + l/2,x 3 )B £i (x° 4 + 1/2,2-4))* x x °>x° ■ 

(4.22) 

The point of the considerations above is that the 
singularities of Gei£ 2 i 3 £ 4 (k) are the same as those of 
DtitztstA&ftk)- 

X 3 ff Xi 



ff = 


Xi 


- x 3 






T = 


X 3 


- X 2 




\ 


f = 


X2 


- X\ 


X\ 


— >- 



Figure 1: Lattice relative coordinates. 

Analogously, reinstating the linear combinations, for 
Mi^i 2 i 3 i 4 , we obtain the spectral representation 

M eiS2S3Si (k)= M SlS2 s 3 sM - (2^) 3 J\ J T3 f(k°, A ) 

x8(k- A) d x od x {? Sl (l/2,ui)f a2 (l/2,u 2 ), 
x£(A°,A)f S3 (l/2,U3)r S4 (l/2,«4))-H, 

where, for x° ^ 0, and with x = (x° = v a — u , x) € Z , 

M slS2S3S4 (k) = ^ ] e M slS2S384 (x =0,a;). 

The singularities of M S1S2S3S4 , for k = (k = = 
0) and e ±x < 1 are points in the mass spectrum. To 
determine M SlS2S3S4 we employ a lattice version of the B- 
S equation, as discussed below. To detect baryon-baryon 
bound states below the lowest two-baryon threshold, we 
write k — i(2fh — e), where ffi is the smallest baryon 
mass and determine the value of e > that gives the 
singularity. We call this special value of e the binding 
energy of the corresponding bound state. 



4.3. The B-S equation, the ladder approximation 
and 1 = bound states 

Next, we write a B-S equation for M as, 

M = M + MqKM , (4.23) 

and in contrast to the continuum case, we use what we 
call an equal time representation. 

In terms of kernels, with 2;^ = x\ and 2:3 = x\, we have 

M(xi,X2,X 3 ,Xi)=M (xi,X 2 ,X 3 ,Xi)+jM (x 1 ,X2,yi,y2) 

x K(yi,y 2) y 3 ,yi)M(y 3 ,yi,x 3 ,Xi) 
x %i ~ 2/2)%° - yl)dyidy2dy 3 dy i , 

where we use a continuum notation for sums over lat- 
tice points. if(yi, 7/2, J/3, 2/4) is called the B-S kernel. We 
remark that in the imaginary time continuum quantum 
field theory formulation with Euclidean invariance all the 
B-S equation coordinates are taken to run over all space- 
time. 

Formally, this B-S equation defines K = Mq" 1 — M~ x . 
Our strategy is to compute K in the leading k order, 
which we call a ladder approximation L to K, and then 
solve the B-S equation for M using L. In the B-S equa- 
tion, Mq is obtained from M by erroneously applying 
Wick's theorem to the one-baryon composite fields B in 
M such as, 

- (B £l (xi)&£ 3 (x 3 )) (Bi 2 {x 2 )B ti (xi)) 
+{B ei (xx)B e4 (xi)) {B e2 (x 2 )B l3 (x 3 )) , 

for the term {Bi 1 (x 1 )Bt 2 {x 2 )B l3 {x 3 )Bi 4 {xi)) in 
Mi 1 i 2 i 3 i 4 (xi, x 2 , x 3 , X4). With this, we obtain 

Fo,e 1 e 2 e 3 e4 ( x i > x 2 , 2:3, x 4 ) = -G^ 3 {x x , x 3 )Gt> 2 i 4 (x 2 , x A ) 

+Ge 1 t 4 (xi,Xi)Gt 2 t 3 (x2,x 3 ) . 

(4.24) 

By referring to the definition of the two-baryon function 
Gi 1 £ 4 of Eq. H3.10fl . we see that Eq. I|4.24|l holds for all 
temporal orderings. 

In terms of the (£, ff, r) lattice relative coordinates and 
taking the Fourier transform in r only, the B-S equation 
becomes (see [2^1 

M(i, ff, k) = Mo(l rf,k) + f M (i, f'> k) 

x K(-g', -ff', k)M(rj', ff, k) di'dff' . 

The above B-S equation corresponds, roughly speak- 
ing, to an equation for the kernel of a one-particle lattice 
Schrodinger operator resolvent equation. The B-S kernel 
K(—£', —ff 1 , k), in general, acts as an energy-dependent, 
non-local potential (see Refs. [UE3). M is like the re- 
solvent of the dressed noninteracting particles and M is 
the interacting resolvent. To facilitate the analysis, it is 
desirable to have K with the shortest range as possible. 
For example, in weakly coupled scalar models, with lo- 
cal interactions and in the weak coupling regime, as well 
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as in other spin and stochastic systems, the dominant 
contribution to K is local, and the ladde r appro ximation 
dominates in the full model (see Refs. E3, In 
general, the hyperplane decoupling method is used to ob- 
tain decay properties of K, which allows us to control the 
solution to the B-S equation. 

To make sense of the B-S equation and the kernel K, 
we now define the space where M and Mq act as matrix 
operators, which is given by 1' 2 = £ 2 ® £ 2 , such that / = 
/i + h G ^2) A C/2) is associated with the first 13 (last 
7) components; /1 € £ 2 (A S ), f 2 6 £ 2 (A a ), where s ( a ) 
means the symmetric (antisymmetric) subspace and 



A s = {( Xl ,x 2 ;j) e Z 4 x Z 4 x (1,2, ...,13)1^ = ^}. 



and 



A a ={( Xl ,x 2 ;j) eZ 4 xZ 4 x (U,...,20)\x° 1 = x° 2 ,x 1 ^x 2 }. 

Using the orthogonality properties of the two-point 
functions, in the total spin basis, we have 



M, 



o°s,Sz,S',S' (xi,x 2 ;x3,X4,) = S S s'5s z s' z [-S(x 1 - x 3 ) x 



S(x 2 - x A ) + (-if x 
6{x x - x 4 )S(x 2 - x 3 )} 



and the same for M^ Tz T , T ,, T = 0,1,2,3. Hence, 

Mq -* acts as —21 7^ on £' 2 , and Mo is invertible. To 
define the B-S operator K = Mq — M~ x , we also need 
M~ l . For M(xi, x 2 , x 3 , x 4 ), setting 

M(0) = M(x, x, x, x) = M (0, 0, 0, 0) , 
we find that the k° order is given by 

(xi, x 2 ,X3,X4)= M^(0)S(xi - x 2 )S(xi - x 3 ) 
x 5{xt - Xi) + M^ 0) [l- 
S(xi - x 2 )5(xi - x 3 )5(xi - xa)] , 

such that M(°) agrees with Mq°\ except at coincident 
points. Hence, to show M is invertible we need the in- 
vertibility of the matrix M^(0) which we analyze next. 

By a lengthy computation given in Appendix C, we 
find that, in the total spin basis, 



m(°> (o) = m[°> (o) ffi m\° ] (o) e m[ 0) (o) e m^ ] (o) 

where 

1 V ' 9 1 -8^5 16 



(4.25) 
(4.26) 



and M^^O) acts as —41 on the components with indices 
from 7 to 13. The eigenvalues of M^°\o) are —4 and 
0. The null space prohibits the invertibility of M^°\0). 
1 ( 2 \ 

The zero eigenvector is -| I I and the first line of the 
eigenvalue equation is 



en 



T T=1 , T J)(°)=0. 



The vanishing of this average suggests that there is a 
pointwise linear dependence relation at coincident points. 
Indeed, this turns out to be the case, and is proven in 
Appendix C that, at coincident points, 

T s =i,s z + ^Tt=\,t z = , (4.27) 

for S Z =T Z = -1,0,1. 

The truncated four-point correlation M still has a null 
space in the symmetric subspace of £' 2 . For this reason, 
we are forced to reduce the symmetric subspace for coin- 
cident points only, and we take the restriction of M<-°\0) 
to the orthogonal complement of the null space in the 
symmetric subspace. So, for coincident points we have 
a 10 dimensional space. For non-coincident points, M^> 
agrees with Mq -* ^ 0, so that K is well-defined. 

With this definition of K, and using its decay proper- 
ties shown in Appendix D, there is a local potential, and 
energy-dependent space range-zero potential, as well as 
an energy-independent range-one potential. As in Ref. 
, we modify M at coincident points by a multiplicative 
constant 1/2. Precisely, we replace M(x±, x 2l x 3 , X4) by 

M'(x ll x 2 ,x 3 ,x i ) = h(x 1 ,x 2 )M(x 1 ,x 2l x 3 ,x 4: )h(x 3 ,X4) , 
where 

h(x, y) = [l- S(x - y)} + ^5(x - y) , 



and the B-S equation becomes 



M' = M Q + MqK'M' . 



(4.28) 



The r variable Fourier transforms M and M' clearly have 
the same singularities. For simplicity of notation, we 
drop the prime and keep using M, M, K and K. 

Next, to define our ladder approximation, we write a 
Neumann expansion representation for K = Mq — M _1 . 
With 



5M = M — M (0) 



SMn = Mn - M, 



(0) 



we have 



K= Y.n^i-miM^SMoYM^- 
-{M^^SMfM^- 1 }. 



(4.29) 



For small k, the dominant contributions to K come 
from coincident points at n = and a space range- 
one potential of order k 2 . As discussed in Appendix 
E, it turns out that, because of cancelations K^°\0) = 

(M^VHO) ~ (M (0) ) _1 (0) = 0. Surprisingly, there is an 
attractive gauge field correlation effect which is respon- 
sible for this result. More explicitly, as a contribution to 
the definition of h(x,y), we must consider the effect of a 
gauge integral over six overlapping bonds with the same 
orientation, which is given by Iq in Eq. I|D6|) . Moreover, 
it is worth mentioning that if these six gauge fields are 
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uncorrelated three by three, meaning that Iq is replaced 
by (1 3 ) 2 , 1 3 being the integral over three overlapping 
bonds with the same orientation [see Eq. l)B.ljl ]. the as- 
sociated attractive contribution to K^(0) becomes re- 
pulsive, leading to an overall repulsive K^°'(0). 

Going further in the k expansion, since 8M is 0(k 2 ) 
and SM is 0(k 3 ), we find that the dominant contribu- 
tion to K is given by the n = 1 term in the Neumann 
expansion of Eq. I|4.29|) . This defines our ladder approx- 
imation to K and is given by 

L = K^k 2 = (M(°))- 1 (5M)( 2 )(M(°))- 1 K 2 . (4.30) 

To obtain an explicit formula for the kernel of L of Eq. 
(|4.30|l . first we consider the contribution to L coming 
from the non-oriented geometric bond linking the sites 
and e . It is given by 

K^(0, e 1 ; 0, e r )+KW(0, e 1 ; e 1 , 0) = (M^)~\0, e 1 ; 0, e 1 ) 
x M( 2 )(0, e 1 ; 0, e^ikfC ))" 1 ^, e 1 ; 0, e 1 ) 
+ (M(°))- 1 (0, e 1 ; e 1 , 0)M(2)(0, e 1 ; e\ 0) 
x (M(°))- 1 (0,e 1 ;e 1 > 0). 

Since M^ ' = = —2 1, for non-coincident points, 

each (M(°)) -1 factor can be replaced by —1/2. In the 
total spin basis, taking into account the space coordinate 
exchange symmetry (anti-symmetry) of T for S = 1 , T — 
1, 3 (S = 0, T = 0, 2) given in Eq. |@2TJ|, the second term 
is the same as the first one. 

The contribution (0, e 1 ; 0, e 1 ) n 2 is obtained by ex- 
panding the hopping term exponential in the action ap- 
pearing in the numerator, which brings down a pair of 
oppositely oriented overlapping gauge field bounds. The 
gauge field integral of Eq. (|D7|I is performed using the 
Peter- Weyl theorem (see Ref. for other more com- 
plicated and useful gauge integrals and Appendix D), 
and gives a factor 1/3. In this way, the 11 element 
MWfO.e^O,^) of the matrix is 

where, in the individual spin basis, 
_ — , 11 11 

' — \ ' r 2 2 r 2 2 V 

^1,83,83,84 — / , C I 1 l2 C hl4 

(Bi IlSl Bi hS3 ^ aiaif ^ aiP2f2 )^ x 

(Bi l2S2 Bi hSi i> a2a2f2 ip a2f3l fl ) (°) T^ ifJi T~f 02 , 

(4-31) 

with all fields at a same site and where the prime here 
means the sum over with Ji + 1% = = 

J3 + 14. A similar expression holds for the compo- 
nents 12 {sxsrtzt^), with (1/2, ^3, S3) and (1/2, J4, S4), 
respectively, replaced by (3/2, J 3 ,i 3 ) and (3/2, J 4 ,i 4 ); 
with (1/2, Ji,s 2 ) and (1/2, 1 2 ,s 2 ), respectively, re- 
placed by (3/2,Ji,ii) and (3/2, J 2 ,t 2 ) for the compo- 
nent 21 (iii 2 s 3 s 4 ); and with (1/2, /i,s x ), (1/2, J 2 ,s 2 ), 
(1/2, /3,s 3 ) and (1/2, 14,34), respectively, replaced by 



(3/2, Ji,h), (3/2, J 2 ,t 2 ), (3/2, J 3 ,t 3 )and (3/2, J 4 , t 4 ) for 
the component 22 (tit 2 t 3 t4). 

Taking into account the other directions and using 
translation and rotation invariance, we obtain for the ker- 
nel of Jf' 2 ), in the total spin basis 

K^(xi,x 2 ;x 3 ,x 4 ) = jgvW J2<r,j S ( x 2 ~ x l -oe 3 ) x 

S(xi — x 3 )8(x 2 — xa)8(x 2 — x 3 ) , 

(4.32) 

where here v^> is the total spin basis matrix. 

By a lengthy computation (see Appendix E), in the 
individual spin basis, the symmetric matrix i/ 2 ) — 
[ v ij]i,j=i,—,20 has the following nonzero elements: v\\ — 

-1/3 = W44; "17 = -4/V6 = "113 = "412 = ^418! "110 = 

4^2/3 = W415; "22 = 10/3 = v 33 ; v 23 = -11/3; v 2 $ = 
2V2 = v 31T ; wa U = -4V2/3 = v 314 ; Wsii - 2^2/3 = 

"2 14! "55 — -3 = "20 205 U66 = -2 = "99 = "16 16 = "19 195 
"69 = -1 = "77 = VgU = U1212 = "13 13 = U1417 = 
"16 19 = "18 185 V 7W = -2\/3 = W1013 = "12 15 = "15 185 

vio 10 = -5/3 = wis 15; "11 11 = -4/3 = vn 14 = V14 14. 

It turns out that has an especially simple, almost 
diagonal form, in the total spin basis. Using the basis 
transformation A given in Section 4.1, it is written as 

v<® = W ®W ®W ®W ld ®W 2 ®W 2d , (4.33) 



with W = - /— - I (with eigenvalues ±3), 

3 \ -vol) 1 J 

W 2 = f J I (with eigenvalues —1 and 9), W\d — 

— 3Pi3, W 2 d = —P22, where P13 is the orthogonal pro- 
jection on the 7 th through 13 th components and P 22 on 
the 16 th through 20 th . 

The interaction present in Eq. I|4.31J) arises from a 
quark-antiquark exchange but, by inspection of the me- 
son particle fields (see Refs. [a,ll3)i this is not a meson 
particle exchange. To see this, consider the field averages 
in Eq. 14.31fl . The spin indices of the B and B fields are 
both lower. This forces the spin indices of the ipij) fields, 
in each of the two averages, to be either both lower or 
both upper. A true meson particle has one upper and 
one lower spin index. We refer to this interaction as a 
quasi-meson exchange. 

Using the total spin basis representation for v^ 2 \ we 
return to the determination of bound states and solve 
the B-S equation using L (see Refs. [jj El for details). 
In terms of relative coordinates, we have the space range- 
one energy-independent potential 

L(H, n, k°) = K V 2 > J2 £>(*e* - , (4.34) 
and the B-S equation, in the ladder approximation, has 
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the solution, 

(4.35) 

where the Mo lattice indices are restricted to the spatial 
nearest neighbors of zero, and we suppress the k° de- 
pendence and spin indices on the r.h.s. Reinstating the 
spin indices, the resulting matrix is 120 x 120, where we 
have implicit matrix multiplication in the spin indices Mq 
means Mo with indices restricted by the multiplication 
of v*- 2 -* and we suppress k° dependence on the r.h.s. 

To go further, we need properties of Mq. In Appendix 
F, using the spectral representation for the two-baryon 
functions in Eq. I|4.24|l . after passing to lattice relative 
coordinates, we derive a general spectral representation 
for Mo- Assuming that the two-baryon function is diago- 
nal in the spin indices in the individual spin basis [which 
holds up to and including 0(k 5 ), see Ref. we show 
that the dominant contribution is given by, in the total 
spin basis, 



M, 



o,SS.S'S>£l ff, fc°) = [-2(27r)- 3 / T3 G(p)G(p)dp 
-2(2tt) 3 £ £ J T3 f(k°, A°A' )d A oa ? (A ) 



x d X 'oap(X'°) S(p, £, rf)dp 5 S s> $s, 



(4.36) 



where 

s (p> = £si cos p4 cos p.rf + 6 so sin p.£ sin p.ff ; 

and the same for T = 1,3 (T = 0, 2) replacing 5 = 1 
(5 = 0). 

To continue, we set k° = i(2m—e), where we recall that 
m denotes the minimum of the baryon masses and e > is 
the bound state binding energy. Recalling the definition 
of p\ given in Eq. I|3.14JI , we further make the following 
approximations: a) we retain only the product of one- 
particle contributions; b) w{p) — m k 3 p 2 /8; c) m ~ 



-31n« and d) Z{p) k, -(2tt) 
The approximate Mq becomes 



-(27T) 



Mo 



SS Z S'S' 



(27T) 



' Jt 3 



1 



1 



1 



x Z(p,£,ff)dp 5 S s> S Sz si ■ 

(4.37) 

With all these approximations, the singularities of 
M , below the two-baryon threshold 2m, occur as ze- 
roes of the determinant of the 120 x 120 matrix 1 — 
K 2 u( 2 )Mo(<Tie 51 , 02 e 5 ' 2 ). Anticipating that the binding 
energy e is of order k 2 , we also approximate the denom- 
inator in Mq to get 



Mi 



o,ss z S 



>sdt,ff,k )> 



(27T) 



(1 



1 ' 3(fi,t,v)dp, 
/t 3 



and the same for T's instead of 5's. These approxi- 
mations are consistent with finding a bound state solu- 
tion with binding energy greater than (D(k 3 ). Note that 
Mo(<7ie J1 , <T2e n ) is diagonal inji andj2, and that the in- 
tegral takes the values ±(27r) 3 /2. Thus, the bound state 
condition becomes 



det[l + 4(l-e" e ) 



-V 2 Vl 



0, 



(4.38) 



for a 20 x 20 matrix. Next, using the spectral representa- 
tion for y( 2 \ if A is a negative eigenvalue, then the bound 
state condition is 

i +4A(i - e - e rv = 0, 

with approximate solution e = 4|A|k 2 , which is the bound 
state binding energy. Hence, we have a bound state for 
each negative eigenvalue of v^- 2 \ with multiplicity given 
by the space dimension 3. For these bound states we 
have attractive potential wells at ie^, j = 1, ...,d = 3, 
which would give us a multiplicity 2d. However, we only 
have symmetric or antisymmetric wave functions, which 
reduces the multiplicity to d = 3. 

Referring to the spectral decomposition of t/ 2 ) in the 
total spin basis, and recalling the ordering of the total 
spin basis, we see that the most strongly bound, bound 
states (A = —3) are associated with a superposition of 
p — n and A — A total spin 1 states; and also the A — A 
total spin 3 states. The more weakly bound, bound states 
(A = —1) are associated with the superposition of p — n 
and A — A total spin states; as well as the A — A total 
spin 2 states. 

We now explain the superposition of two-baryon states 
occurring in the bound state associated with the super- 
position of p-n and A — A. Similar considerations apply 
to the other bound states. Going back to the decompo- 
sition of given in Eq. (|4.33ll . this state is associated 
with the negative eigenvalue —3 of one of the three 2x2 
W's. Consider the first W and label the basis elements 
(5,5 Z ) = (1,1), (T,T 2 ) = (1,1) by = 1,2. The corre- 
sponding eigenvector is W-3 — (v5/3, 2/3). If we start 
from the B-S equation for (io_3, Mw-3), it reduces to a 
single equation with i/ 2 -* substituted by —3. Here, we 
have used the fact that M commutes with in our 
approximation. But (w_ 3 ,Afw_ 3 ) can be written, for 



,(0) _ ^(o) 



-«J) _ -(0) 
L 3 — x 4 I 



as the average 



( [(V5/3)T^l St=1 ( Xl ,x 2 ) + (2/3)T*Ll T ^ 1 (x 1 ,x 2 )} 
x l(V5/3)T^l s Jx 3 ,x 4 ) + (2/3)75 



T=l,T i = l( 2 '3, Xi)] ), 



SO 



that the superposition of two-baryon states 
given by the field [(V5/3)Tgi 2 1Sz=1 (x 1 ,x 2 ) + 

(2/3)7^.^ T =1(^1,^2)] more closely describes the 
field which creates the bound state. 
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4.4. Bound states for the I — 3 sector 

In this subsection, we show that the application of our 
method to the maximum total isospin 1 — 3 sector (freez- 
ing all isospins to +) leads to no baryon-baryon bound 
states if S = 2, 3, but bound states do appear if S = 0, 1. 
The analysis here parallels that of the I — case, but it 
is algebraically much simpler, and we will be brief. 

The one-baryon states are associated with the fields 
B t = Bs3 t , t = ±3/2, ±1/2 of Eq. j£5jL The two- 
baryon states that we consider here, in the total 1 = 3, 
I z = 3 sector, are those generated by the product of 
the above one-baryon fields f tlt2 (xi, x^), x\ = x^. The 
bound states are formed with a quasi-meson exchange k 2 
space range-one potential which dominates and, in our 
ladder approximation, the space range-zero potential is 
zero. 

After interchanging spin and isospin, the calculation of 
this potential is the same as for the lattice QCD model 
in 2 + 1 dimensions and with 2x2 spin matrices and two 
flavors (see Ref. 0), and the bound state spectral results 
are the same, apart from the additional space dimension. 

In the total spin basis, the two-baryon states are given 
by 

T TTt (xi,X2)= Yl ^^1*2(^1^2), 

ti+t 2 =T z 

where c^T* are the C-G coefficients for the coupling of 
two spin 3/2 states to give total spin T = 3,2, 1,0, and 
z-component T z of total spin. They obey the symmetry 
property 

T TTz (xi,x 2 ) = (~l) T T TTz (x 2 ,xi) , 

for T = 3,1; for T = 0, 2 (1,3) they are space sym- 
metric (antisymmetric). So, at coincident points, they 
vanish for T = 3,1; for T = 0, 2 (1,3) they are sym- 
metric (anti-symmetric). Also, from the pointwise rela- 
tions BisB+s = BijBii = BisLi = B. iB . 1 = 0, 
- 1 - 2 2 2 2 2 ' 2 2 2 

we find that Tt =2 t z is also zero at coincident points. 
A computation gives the coincident point linear depen- 
dency relation Bs_B_3 = —BiB_x, and we find that 

22 22 

Tqq = 2BsB_3. 

The four-point correlation in the unmodified, shifted 
form, and in the individual spin basis, is taken as, with 
x\ = x%, £3 = £4, and for x\ < x 3 , 

Mt 1 t 2 t 3 u(x 1 ,X2,x 3 , £4) = (B tl (xi)B t2 (x 2 )B t3 {x 3 )B ti (x 4 )) , 

and with the usual interchange of barred and unbarred 
fields and taking the complex conjugate of the average, 
for X® > x 3 . The associated uncorrelated or Wickified 
four baryon correlation is, for x\ < x®, 

Mo ttl t 2 t a u(xi, x 2 , x 3 , x A ) = 

- (B t ,( Xl )B t3 (x 3 )) (B t2 {x 2 )B tl {x 4 )) 
+ (B tl { Xl )B u (x 4 )) (B t2 (x 2 )B t3 (x 3 )) . 



The four-baryon correlations in the total spin basis are 
obtained by taking the corresponding C-G linear com- 
binations. For T = T = T z = T' z = 0, we obtain 

M^ T , T ,(0) = -4 and M ( °^ T , T , (0) = -2. 

As in Section 4.3, an /i-modified B-S equation 
and operator kernels are introduced. Then, the 
energy-independent space range-zero potential K^°> = 
[M^]- 1 - [AfW]- 1 = 0. Analogous to Eq. jOUJl , the 
dominant interaction is a k 2 quasi-meson space range- 
one potential which we take as our ladder approximation 
L = K^n 2 to K. Here, with and v {2) in the total 
spin basis, as in Eq. (|4.32|l 

K^(x 1 ,x 2 ;x 3 ,x 4 ) = jzV i2) J2a,j S ( x 2 - xi -ae j ) x 
S(xi — x 3 )S(x 2 — Xi)S{x% — x 3 ) , 

where, in the individual spin basis, 

V iut 2 ,t 3 M = r ai/3i r a2 e /3 2 ( B tl B t 3 ^a iaifl lpa lf 32f2) {a) 
X (S t2 i3 t4 ^ Q2Q2/2 ?/'a 2/ 3 1 / 1 ) (0) ■ 

(4.39) 

Next, using similar arguments to those in Appendix E, 
we can write 

v tut 2 ,t 3 ,t4 = ( B ti B t 3 : ^axccfAaifSh : > (0) 

where in the expansion of the k = average using Wick's 
theorem the Wick ordering : : forbids the contraction be- 
tween the explicit internal if) and tp fields. This restric- 
tion implies that fi = f 2 = +■ In order to calculate 

(2) 

v ti ' t t t , we also take into account spin constraints and 
use conjugation and the k — spin flip symmetry. 

Doing this and suppressing the + isospin indices of the 
ip, it is enough to calculate 

(B1B1 : $ a +ipa+ :) (0) = -2 = 2(BiBi : ^ a -^a- :) (0) , 
(BzBi :$„+1> a - :)(°) =-V3, 
(BsBz :$ a+ i; a+ :)(°) =-3, 

and : ^ a+ Va- :) (0) = -2. 

In this way, in the individual spin basis, using the 
same ordering as in the zero isospin case, for the 16 x 16 
submatrix, the symmetric matrix = [%']ij=i „,i6 
has the following nonzero elements: vn = 9 = Wi6i6! 
v 22 = 6 = W55 = W1212 = U1515; v 25 = 3 = v 33 = v i7 = 

V88 = «99 = V1013 = U1215 = VU14,; V 3 Q = 2V3 = t<69 = 

van = vnu; V66 — 5 = uim; v 77 — 4 = v 7w = v W w- 
The matrix is easily diagonalized in the following re- 
ordered basis: 3, 6, 9, 4, 7, 10, 13, 8, 11, 14, 2, 5, 
12, 15, 1, 16. Letting V denote this matrix, we have 
V = Vi © V 2 ® Vi 9 V 3 9 V 3 e V 4 e V 4 , with Vi,2,3,4 sym- 
metric, and where Vi is 3 x 3 with nonvanishing elements 
Vi-u = 3 = Vi;33, Vi ; 22 = 5, Vi ; i2 = 2\/3 = Vi ;23 ; V 2 
is 4 x 4 with nonvanishing elements V2 ; 22 = 4 = V 2; 3 3 = 
V2;23, V 2 ;i2 = 3 = V2 ;3 45 V 3 is 2 x 2 with nonvanishing 
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elements V3 ; ii = 6 = ^3,12 — 3; V4 is 1 x 1 with 

nonvanishing elements V^-n = 9. 

Concerning their eigenvalues and eigenvectors (T, T z ) 
(the total and z-component of the total spin), we have: 
9 (3, 1), 3 (2, 1) and -1 (1, 1) for Vi in the first spot; 9 
(3, 0), 3 (2, 0), -1 (1, 0) and -3 (0, 0) for V 2 ; 9 (3, -1), 3 
(2, -1) and -1 (1,-1) for Vi in the second spot; 9 (3, 2) 
and 3 (2, 2) for V 3 first spot; 9 (3, -2) and 3 (2, -2) for 
V3 in the second spot; 9 (3,3) for V4 in the first spot; 9 
(3, —3) for V4 in the second spot. 

As we see, in the total spin basis, is diagonal. Al- 
ternatively, it can be also diagonalized by the orthogonal 
transformation relating the total and individual spin ba- 
sis. Similar to what occurs for the 1 = sector, we 
shall see that the negative eigenvalues are associated 
with bound states. Referring to the above, the nega- 
tive eigenvalue —3 (—1) is associated with T, T z = (0,0) 
(T,T Z = (1,1), (1,0), (1-1)). 

Proceeding as in the 1 = case, in terms of lattice rela- 
tive coordinates, the space range-one energy-independent 
potential in the ladder approximation is, as in Eq. I|4.34|) . 
given by 

L{1 ff, k°) = k 2 v {2) £ S(ff- 05(aeJ - f) , 

with v^ 2 ' in the total spin basis and the B-S equation 
has the solution given by Eq. I|4.35|l . Here M (^,f},k°) 
is derived as in Appendix F and in the total spin repre- 
sentation is given as in Eq. I|4.36[) with S, S z replaced 
by T, T z , as well as their primes and with now 

-{p, f , ff) = (S T2 + 6to) cosp.£ cosp.rf 
+ (St3 + Sti) smp.^smp.ff . 

Finally, approximating Mq as before, for each negative 
eigenvalue A of we have the bound state condition 
1 + 4A(1 — e~ e )~ 1 K 2 = 0, with the approximate solution 
of binding energy e = 4|A|k 2 . Hence, we only have bound 
states of spin one and zero, with multiplicities 3 and 1, 
respectively. 



5. CONNECTION BETWEEN THE B-S 
EQUATION AND A LATTICE SCHRODINGER 
EQUATION 

In this section, to understand our spectral results in 
a simple way, we exploit the similarity between our 
lattice B-S equation for M(£, ff, k°) and a one-particle 
Schrodinger resolvent operator equation for the lattice 
Hamiltonian 

H = H + V , 

with a kinetic energy operator Hq and a potential energy 
operator V. 



For k and e small, we further approximate Mo(£, ff, k ) 
of Eq. fl4"37|) by 

2 f 1 

M ° = ~7T^3 / 3 21A I Z(p,t,ff)dp, 

for pi given in Eq. H3.14[) and 2 depends on the total 
spin. 

Consider the Schrodinger resolvent equation, 

{H - z)- 1 = (H - z)- 1 - {H - z)- 1 V(H- z)' 1 , 

for z not in the spectrum of Hq and H. 

Upon comparing with minus the B-S equation in our 
ladder approximation, i.e. —A/ = —M0 — M0LM, we can 
make the identifications 

3 

H = ^(-A e ) , L = V , z = - e -{z<0), 

where Ai is the lattice Laplacian with Fourier transform 
— pj. The Hamiltonian H = Hq + V acts on the space 
£20 (Z 3 ), which decomposes into a direct sum of space- 
symmetric and antisymmetric subspaces depending on 
the total spin indices (omitted here for simplicity of no- 
tation!). For / e l\o{l?), for each x, fn{x) is the Rth 
component where R runs over the 20 total spin basis in- 
dices. The subspace is symmetric (antisymmetric) for 
£,ff) = cosp.£ cosp.ff (sinp.£ sin fi-ff), corresponding 
to total spin S= 1, T= 1,3 (5 = 0, T= 0,2). 
The kernel of V is 

V& ff) = k\W 5(£- aei)5(£- ff) , 

and, for small k, V dominates Hq. 

In this case, the negative energy bound state eigen- 
functions of H with eigenvalue A = z < arc those of V 
which are easily seen to have components 

f { B k \x) = a R [S(x~e k )±S(x + e k )}, 

where a is an eigenvector of the eigenvalue equation 

k 2 v^a = \a . 

The plus (minus) sign occurs for the symmetric (anti- 
symmetric) subspace of £20 (^ 3 )- The eigenvectors are 
localized at ±e fc , respectively. When the free Hamilto- 
nian Hq is not neglected, the negative eigenvectors of H 
get spread out. 

We now reinstate Hq and obtain the spectral repre- 
sentation of H where we write V(£,ff) = 2K 2 fiJ2j <K£ ~ 

e J )<5(£ — ff), where fj, is a negative eigenvalue of v^ 2 \ and 
we are in the space symmetric or antisymmetric subspace. 
As H acts in the space of symmetric or antisymmetric 
subspace of £20 (^ 3 ) the 6(£ + e 3 ) term in the potential 
can be dropped and the term <5(£ — e J ) gets multiplied by 
2. The space symmetry or antisymmetry is incorporated 
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automatically in the definition of (Ho — z) 1 through the 
function £, ff). 

The spectral results below hold for all k > while 
the correspondence with the B-S equation is valid only 
for k <C 1. With the above observations, the resolvent 
equation becomes 

(H - z)- 1 ^, ff) = (H - z)-\l ff) - {H - z)- 1 ^, e>) 
x 2k 2 /j,(H - z)- 1 (e j ,ff), 

with solution 

(H - z)-\lti) = {Ho - z)-\£,rf) - {Ho - s)" 1 ^) 
x 2 K 2 t iB jk (H -z)- 1 (e k ,ff), 

where B = and A = [1 + (H — z)~ x 2n 2 fi] restricted 
to ^({e 1 , e 2 , e 3 }) so it is a 3 x 3 matrix. The matrix A 
has the structure 



A = 



where a = 1 + 2k 2 {i(H - z) _1 (r,t), (r,t) = (e^e 1 ), 
(e 2 ,e 2 ) or (e 3 ,e 3 ), and & = 2k 2 /j,(H - z)~ 1 (r, i), (r,i) = 
(e^e 2 ), (e^e 3 ) or (e 2 ,e 3 ). Furthermore, detA = (a — 
&) 2 (a+2&), and the eigenvalues of A are Ai = A 2 = (a— b), 
A3 = a + 26, with corresponding normalized eigenvectors 
h = (n,r 2 ,r 3 ), f 2 = (h,t 2 ,t 3 ), f 3 = ^=(1,1,1), where 

£ i r i = £ i t i =0and/ 1 l_/ 2 . 

The matrix i? has the same structure as A and is given 

by 




B 



1 



(o-6)(a + 26) 



a + fe -b -b 
-b a + b -b 
-b -b a + b 



The negative energy bound states occur as zeroes of the 
denominator of B, and are given by the zeroes of (a — b) 
and (a + 2b). 

Denoting minus the binding energies by z — —e, for 
d = £2, £3, we have the conditions, for z = —e\ = —e 2 

= 1 +2K 2 f i[(H - z)-V, e 1 ) - {Ho - z)-\e\ e 2 )] , 

and, for z = —£3 

= 1 +2k 2 h[(Ho - z)-\e\ e 1 ) + 2(tf - z)- 1 (e 1 , e 2 )] . 

The corresponding numerators of the resolvents {H — 
z)^ 1 , in Fourier space, are (cosp 1 ) 2 — cosp 1 cosp 2 
and (cosp 1 ) 2 + 2 cos p 1 cos p 2 , in addition to the factor 
E{p,£,ff). Expanding the denominator (p 2 — z) shows 
that £3 > £1. 

Recalling that the resolvent contains all spectral infor- 
mation, we now obtain the bound state eigenfunctions. 
To this end, from the spectral theorem 

hm (-£, - z)(H - z)- 1 ^, ff) = £7=1 1>n@$iM 

Z — > — € i J 

= -{Ho + ar^le^F^Ho + e i )- 1 (e k ,ff) , 



where F^ k — lim z ^_ ei 2n 2 [iBjk and ipij(0 are eigen- 
functions of H with eigenvalue — ef, j = 1, 2, . . . , m t la- 
bels the TO, degenerate states. 
For i = 3, a = —2b and 



1 



1 1 1 



F 3 = — 1 1 1 
3c3 \ 1 1 1 

where, upon expanding (a + 2b) in the denominator of B 
in powers of (z + e 3 ), we have c 3 = 2[{H — z)~ 2 (e 1 , e 1 ) + 
2{H — z) _1 (e 1 , e 2 )]| z= _ £3 . The eigenvalues of the ma- 
trix are and 3, with multiplicity two and one, respec- 
tively. The corresponding eigenvectors are hi, h 2 and 
hs = (1, 1, 1)/V3; in hi and h 2 the sum of the compo- 
nents is zero and hi is orthogonal to h 2 . Noting that 
h 3 jh 3k is the jk matrix element of the orthogonal pro- 
jection on the multiplicity one eigenspace, we have 



^3 (£) V>3 (V) =~ (Ho + ei)~% e J ) h 3j h 3k (Ho + ei )~\e k , ff) , 
C3 

so that the normalized eigenfunction is 

M£) = ^(Ho + ei)- 1 (£,ei)h 3j . 

Similarly, for the doubly degenerate eigenvalue £1 = £2, 
we have 



Fi = 



and the matrix has eigenvalues (multiplicity one) and 
3 (multiplicity two) and eigenvectors h 3 = (1, 1, l)/\/3; 
and hi, h 2 where hi is orthogonal to h 2 and the sum of 
components is zero. Hence, 



£m=i,2^im(£)V'im(??) = — (#0 + £1) 1 (£, e ] ) [hij hi k 

Cl 

+h 2l h 2k }(Ho + e 1 )- 1 (e k ,rJ), 

so that the normalized eigenfunctions -0n and V12 are 
given by 

1 




and 



Mt) = -}=(H + e 1 )- 1 {Z,e i )h lj 



^i2(0 = -7=(ffo + ei)- 1 K,e J ')H-. 
v c i 



which are orthogonal to each other since 

X>n(^i2(D = -^(^o + ei)- 2 (e J ,e fe )/i 2fe , 



and from the structure of (i?o + £ i) 1 > ^1 and /12 are also 
its orthogonal eigenvectors. 
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The eigenf unctions ^11,12,3 (£) have exponential decay 
by the Paley- Wiener theorem (see Ref. |3(j), and in fact 

the integrals of the form (Hq + e) _1 (£, ffj occurring in the 
wave functions as well as in the bound state condition can 
be evaluated expli citly in terms of elliptic functions (see 
Refs. jHIH IS Efl El E2 and also Refs. 
We remark that another approach to the eigenfunctions 
is group theoretical noting that the potential is invariant 
under the symmetric group S3 of permutations of the co- 
ordinate components. We see that the eigenfunction ip^ 
is a basis for the one-dimensional identity representation 
of S3; ipn, 1P12 form a basis for the two-dimensional rep- 
resentation of S3. 



6. CONCLUSIONS AND FINAL REMARKS 

We point out that we can determine the baryon-baryon 
bound state spectrum using the more traditional unmod- 
ified B-S equation l|4.23|) . The solution of the B-S equa- 
tion is more complicated as it deals with a three-part 
energy-dependent potential, as explained in Appendix D. 
Using this approach, we do obtain the same bound state 
results. 

To conclude, if we consider contributions to the B-S 
kernel comprised of linear chains of quark, anti-quark 
pairs they result in an exponentially decreasing potential 
with decay rate —2 In clS for the Yukawa theory. Surely, 
it would be interesting to look for the expected distance -1 
Ornstein-Zernicke like correction to this potential and 
to determine the spin and isospin dependencies of the 
binding energy, and the dependence on the number of 
flavors. Especially, the analysis of the SU(3) flavor case 
can shed some light in understanding bound states when 
strangeness and, consequently, A particles are present. 
Also, there is the problem of determining bound states 
of more than two baryons. However, lattice effects are 
expected to be relevant and unrealistic in determining 
the resulting geometric spatial configuration of possible 
bound states. 

We make some remarks on pitfalls that can occur 
when trying to determine bound states by simulation 
methods. For example, consider the total isospin I = 
3, I z = 3 sector treated in Section 4.4. We have 
found bound states with binding energies proportional 
to 3k 2 and k 2 which are associated with the total spin 
(T,T Z ) = (0,0) and (T,T Z ) = (1, (-1,0, 1)), respectively 
If we look at the correlation (7tt z { x i, X2)Ttt z (^3, £4)) 
with x\ < X3, x\ = X2 and X3 = £4, it is zero for 
(T,T Z ) = (1, (— 1, 0, 1)), as TttX x Xi x 2) is space anti- 
symmetric. Thus, these bound states would be missed! 
Even if we take x\ ^ x*i and X3 ^ £4, the decay rate in 
X3 — x\ would not be 2m — an 2 , a > 0, but the smaller 
one of 2m — 3an 2 since we expect there is overlap with 
the (T,T Z ) — (0,0) state. This overlap occurs as there 
is no spin invariance and associated decomposition into 
disjoint spin sectors, at n =/= 0. Thus it seems to be very 
difficult to see all bound states with small spectral sepa- 



rations. 

Finally, and more importantly, we would like to know 
how the baryon spectrum and the bound state binding 
energies behave near the scaling limit. 

This work was supported by CNPq and FAPESP. Al- 
gebraic computations were checked with Maple and Mat- 
lab, for which we thank Louis O'Carroll. We also thank 
P.H.R. dos Anjos for help in unveiling the time reflection 
symmetry of Appendix B. 

APPENDIX A: Isospin orthogonality relations and 
correlation identities 

Here we obtain isospin orthogonality relations for the two 
and four-baryon correlations used throughout this paper, 
and which hold for all k. These relations also hold for 
spin at K = 0, since the action has a SU(2) symmetry in 
the lower components of the spin index, and spin is also 
conserved under this condition. 

We consider, typically, the two-point correlation 
{4>iimi{xi)4>i 2 m 2 { x 2)) and the four-point correlation 
(4>i imi (xi)(f>i 2m2 (x2)(f>i 3m3 {x3)(f>i i m i (x4)), where we sup- 
press the spin index and display only the I and I z labels. 

Considering first the two-point correlation and sup- 
pressing the gauge and spin indices and site arguments, 
the two-point correlation is given by linear combinations 
in the gauge and spin indices, and linear C-G combina- 
tions in the isospin indices of the basic correlation 

(i>i 1 i>i 2 tl>i a $U$hi>ie) = T 123:456 ■ 

We start from the identity, with arbitrary U £ U(2), 

= ((U 1p)t 1 (U l/') l2 (t / '0) l 3^4^5^6) . 

where the r.h.s. is obtained from the l.h.s. using the 
two-flavor isospin symmetry. We schematically rewrite 
this relation as 

(Tt/ 3 )l23;456 = (t/3T)i23;456, (A.l) 

where we set U3 = U <£> [/ <£> E/ . 

We obtain orthogonality relations for the two-baryon 
correlation by relating this identity to the usual quan- 
tum mechanical sum of angular momentum relations. We 
write U — e l9cr / 2 , / = a/2, where 9a = 9 x a x , y a y , 9 z a z , 
or d x ,y,z £ K. Expanding in 9, we have U = 1 + i9I — 

+ . . . such that 

U 3 = 1 <g) 1 ® 1 + i9(I ®l®l + l(g)J(g)l + l(g>l<g)I) 
~ {I 2 ®\®\ + \®I 2 ®\ + \®\®\®l 2 
+21 ® I ® 1 + 2 i" <g> 1 <g> / + 2 1 ® / ® 7) . . . 

Multiplying Eq. 1|A.1[) by W123 on the left and W123 on 
the right, and equating the coefficients of i9 z and — (9 2 + 
9 2 + 9l)/2, we obtain the identities 

M>123(7 1 7 z )i2 3 ;456V456 = Wi23(/ z T)i2 3; 456V456 i ( A - 3 ) 
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and 

?5l23(T/ )l23;456«456 = Wi2^{I 2 T)i23;456V456 , 

where, with an abuse of notation, I z = I z ® 1 ® 1 + 1 ® 
7 z (g)l + l®l(g)/ 2 and I 2 = (Ji + I 2 + ^s) 2 , and /}, as 
usual, acts on the jth spot of the tensor product. 

We will see that I z (J 2 ) has the interpretation of the 
z— component (square) of total isospin. We can use the 
usual complex Hilbert space notation ( , ) ( as (Tx.y,z are 
self-adjoint) and write the above as 

{w,TI z v) = (I z w,Tv) , (w, TI 2 v) = {I 2 w, Tv) . (A.4) 

Thus, if v and w correspond to eigenfunctions with dis- 
tinct eigenvalues we have orthogonality. To be more 
explicit, using the usual quantum mechanics notation 
for a (f3) for spin up (down) state, and suppressing 
the tensor product notation, an isospin state of total 
isospin | and z— component | (|) are represented as 
the C-G linear combination V33 = acta and vsi = 

2 2 2 2 

-^(aaP + a(3a + /3aa), respectively. 

Taking these states into Eqs. I|A.4|) . and using the 
usual quantum mechanical relations gives the orthogonal- 
ity relations. Furthermore, taking the appropriate gauge 
and spin linear combinations to obtain </>^ mf gives the 
orthogonality relations for the two-baryon function. 

Similar to the I z identity of Eq. (|A.3(1 . we also have, 
with I± = I x ± il y , 

(w,TI ± v) = (I T w,Tv) . (A.5) 

Next, we take v — xim and v — Xi(m-i), where \i£ is 
the usual normalized eigenfunction of total isospin / and 
z— component £. Hence, \ie satisfies I±xn = c t Xi(e±i), 
where 

4 = [i(i + i)-e(£±i)} 1/2 . 

Substituting in Eq. I)A.5|) , with the lower sign, and noting 
that cj = cf_ lt gives 

(xi(£-i),Txi(e-i)) = (xu,Txu), 
which implies the identity 

((/>i(e-i)(x)(/>i(e-i)(y)) = (<t>u(%)<i>u(y)) • (A.6) 

Of course, Eq. I)A.6(I implies that all associated one- 
baryon spectral properties are the same, for / fixed, and 
all I z . 

Let us now turn to the four-baryon functions. The 
basic four-point function is 

(V>n ■ ..1pi 6 i>i 7 • ••V'tia) = T X .. .6:7. ..12 ■ 

As before, we have the identity 

TU 6 = U 6 T, (A.7) 



where = U. Decomposing XJq as U3 <8 U3 and 

using similar considerations as above, we obtain orthog- 
onality relations for the z— component and square of the 
total isospin by the appropriate C — G linear combina- 
tions of the product of the one-particle isospin states 4>i£- 
Concerning the Wickified four baryon function Mq of 
Section 4.3, the basic correlation is 

T ,l... 6 ;7...12= -(ipl ■ ■ -^3^7 • ■ ■V'sOv^t ■ ■ -^6^10 ■ ■ ■ ^12) 
+ (lf) 1 . . . -03^10 ■ ■ ■ ^12)<^4 ■ ■ ■ *l>6$7 • ■ ■ V>9>- 

We have the identity 

-fa...^$U) 7 ...{$V) 9 ) (^••■V'e 

x $u) 10 . . .$u) 12 }+ . . .Mi>u) w . . . $U)i 2 ) 

x <v> 4 . . . Mi>u) 7 . . . $u) 9 ) = -{{Ui>)i . . . (c/v) 3 

X ^7 ... l/ig ) {{Utp)i . . . {Ulp) 6 tp w . . . -0i 2 ) 

+{(Ui/>) 1 . . . (U^w ■ ■ ■ ^12) ((Ui>) 4 . . .(Ety)o 
X i>7 . . • ^9) , 

which can be compactly written as 

T U 6 - U 6 T . (A.8) 

We conclude by remarking that the transformation 
properties for T and T of Eqs HA.7|) and (|A.8fl are the 
same. Consequently, we have the same orthogonality re- 
lations for the z— component and square of total isospin 
of two-baryon states. 

APPENDIX B: The hyperplane decoupling method 

In this Appendix, we use the hyperplane decoupling 
method to see how the twenty baryon particles arise 
in the energy-momentum spectrum for our SU(3), two- 
flavor lattice QCD model. It is understood here that 
we first work in a finite-volume model and then with 
bounds on the normalized correlations which are uniform 
in the volume, obtained via a polymer expansion, the ar- 
guments can be extended to the infinite volume theory 
in a standard way (see Ref. 0, |U E3)- This step is 
not treated here, and we work directly in the infinite vol- 
ume. Also, we explicitly discuss the hyperplane decou- 
pling method applied to the time direction, the treatment 
for the three space directions being similar. 

In order to give a brief description of the method, we 
recall that the starting point for the time hyperplane de- 
coupling method is to consider the model with interpolat- 
ing parameters replacing the hopping terms giving cou- 
plings in the temporal direction. In terms of the model 
action, for the hopping terms connecting the hyperplanes 
x° = p and a; = p + 1, the hopping parameter k is re- 
placed by a real parameter n p . (Note that setting k p = 
erases all connections in the time direction between p 
and p + 1 !) (see Ref. 46]). Doing so, as in Sections 3 
and 4, we define K p -dependent correlations e.g. the k p 
dependent zero, two and four-baryon correlations, which 
provide a continuation to complex k p afterwards. 
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The main point here is that decay properties of these 
correlations in direct lattice spacetime are related to 
the vanishing of the first n derivatives with respect to 
k p , at k p — 0. In our case, the first two temporal 
hyperplane derivatives of G vanish, which leads to a 
— 31n/c decay rate. The third derivative has an im- 
portant product structure which implies that the first 
four derivatives of the G convolution inverse T are zero. 
This in turn shows that T has a faster decay rate of at 
least — 41n>t. Associated to these properties, by Fourier 
transform, there are analyticity properties in energy- 
momentum Fourier dual variables. The Fourier trans- 
form r(p) is analytic in a larger p° strip than G(p) so 
that f- x (p) = cof*[f](p)/detf(p) provides a meromor- 
phic extension of G(p) to the larger strip. As seen in 
Section |3J the dispersion curves are determined as the 
zeroes of detr(p° = iw(p),p), so that they are isolated. 
This is a brief account of how, using the spectral repre- 
sentations for the Fourier transformed correlations, the 
method leads to spectral properties. 

Next, we show how the 20 baryon particles arise in the 
two-flavor case, using the hyperplane decoupling method. 
Also, we show how to classify these states according to 
total isospin I, z-component isospin I z and spin s z . Let 
e W a f= £ a bJi>aaif 1 i>ba2 fai>caa f s , which is invariant under 
the exchange ctjfj <-> fk&k, which we call the totally 
symmetric property. As seen below, we can restrict our 
analysis to the lower spin components, ctj = 3, 4. 

For x° < y°, we consider a generic two-baryon correla- 
tion 

G- r r (at, y) = (e^l ? (x) e^l r (y)) . 

ai/ia2J2 V ' aJ X r 5i/i V ' ^ a^h 

We now consider the derivatives with respect to the 
hyperplane parameters k p , at k p = 0. It is easy to see 

that G^Ppl'ty = because of imbalance of fermion fields 

"1/102/2 

ip. For the third derivative, after performing the gauge 
integral using [see Ref. @] 

^3 = JU(g) aibl U(g) a2b2 U(g) a3b3 dfx(g) 

= J 9a 1 b 1 ga 2 b 2 9a 3 b 3 dn{g) (B.l) 
_ 1 , , 

g t aia2a a ^bib2b 3 > 

we get, with w G Z 3 , 

x (e^lAp+l,w)e$l -(y)) (0) , 

(B.2) 

where the structure of r +e restricts the spin indices 
in /? to be lower. Hence, the product structure of Eq. 
(IB. 2(1 is obtained if we take d?i and 6li with only lower 
indices. This formula is also closed meaning that the 
correlation on the l.h.s. of Eq. (|B.2|) is the same as the 
two correlations appearing on the r.h.s. It is important 
to stress that, if we start with a gauge invariant aver- 
age with generic fields with an odd number of fermions, 



(H(x)L(y)) , we do not need to make an a priori guess 
of what the baryon fields are. They emerge as fields to 
achieve the product structure, and is consistent with con- 
finement as the baryon fields are gauge invariant local 
fields with three quarks. 

As explained above, the product property in Eq. (|B.2() 
leads to the existence of particles with isolated dispersion 
curves, after eliminating the redundancies due to linear 
dependencies. 

As each a / takes four values, there are 64 pos- 
sible a f's but, by identifying states which are re- 
lated by the totally symmetric property, this num- 
ber is reduced to 20. The list of 12 representa- 
tives out of the 20 states with their multiplicities in 
the intermediate sums of Eq. (|B.2|) is the following: 
/ = (+,+,+),(-,-,-), a = (+,+,+),(-,-,-), 
multiplicity 1; / = ((+,+,+),_(-,-,-)), a = 
(+,+,-), (-,-,+), multiplicity 3; / = (+,+,-), d = 
(+,+,+), (-,-,-), multiplicity 3; / = (-,-,+), a = 
(+,+,+), (—,—,—), multiplicity 3. The other 8 states 
are given below. 

Eight of the above 12 states with / = 
(+,+,+), (-,-,-) have I = 3/2 and I z = 3/2, -3/2, 
respectively, and s = ±3/2, ±1/2. For the remaining 
four states, we use the fact (from properties of the 
C-G coefficients) that / = 3/2, I z = ±1/2 states 
formed from tensor products of three isospin 1/2 
states have equal weights and are symmetric under 
interchange of any of the individual isospins. In this 
way, these states with d = (+,+,+), (—,—,—) have 
/ = 3/2, I z = ±1/2. Thus, we are left with 8 states, 
with \I Z \ — 1/2 and \s\ = 1/2, given by (including the 
respective multiplicities on the right) 



1)/ = 


(+,-+) 


, d= (- 


h-,+), Iz = 


I a — I 3- 

2 ' 2 ' ' 


2) f = 


(-+,+) 


, d = (- 


h = 


i 8= 1 6- 

2 ' 2 ' ' 


3)/ = 


(+,-,+) 


, d = (- 


K— ,-), h = 


i 8= -i 6" 

2' 2' U ' 


4)/ = 


(-,+,+) 


> d = (- 


h- -),4 = 




5)/ = 


(+,--) 


, d = (- 


h = 


2 ' 2 ' u ' 


6)/ = 


(-+,-) 


, d = (- 


K-,+),4 = 


2 i * — 2 ' ' 


7)/ = 


(-,+,-) 


, d= (- 


h+,-), /, = 


_I a - -I •? 

2 ' 2 ' 


8)/ = 


(+,--) 


, d= (- 


",+,-), Iz = 


2 , 2 , u 



We identify linear combinations of these 8 states which 
we show to have 7 = 3/2 and I = 1/2. The proton 
(neutron) is identified with the I = 1/2, I z = 1/2 (-1/2) 
state. 

Coupling the first two isospins of 1) and 2) (by an- 
tisymmetric linear combination), to give a zero isospin, 
than coupling the third isospin gives the p + state. Sim- 
ilarly, antisymmetric linear combination of 3) and 4), 
gives p_. n + (n_) is obtained by antisymmetric linear 
combination of 5) and 6) ( 7) and 8) ). For the I = 3/2, 
I z = 1/2, s = 1/2 state with a — (+,—,+) take the 
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linear combination 1) + 2) + 2) and note that the second 
2) is not changed by a±fi <-*■ f^a^. Thus, the resulting 
state is totally symmetric in any ctjfj «-> fk&ki hence 
it is a / = 3/2 state. Similarly, the remaining three 
states are obtained with (I,I z ,s) — (3/2,1/2,-1/2), 
(3/2,-1/2,1/2), (3/2,-1/2,-1/2), and are identified 
with the linear combinations 3) + 3) + 4), 5) + 5) + 6) 
and 7) + 8) + 8), respectively. 

The above 20 states are unnormalized ones. Using 
the twenty normalized states of Section 3.1, and their 
multiplicities, and rewriting the intermediate states of 
Eq. I|B.2() in terms of normalized states, we obtain, for 
x° < y° and w € Z 3 , 

(B„,.(s)S«,.'(y)> (3) = -£<iW*)^M> (0) 

w,r 

x(B IIzr (p+l,w)B IIzS ,(y)) {a) - 

(B.3) 

Similar considerations apply for x° > y , and we fi- 
nally obtain the product structure for the two-baryon 
correlation of Eq. (|3.1U|) 

g?X v) = - E G *M X > (p. *)) G «i + y) • 

(B.4) 

Furthermore, for k = 0, using Wick's theorem in the ^ 
fields, we obtain the normalization property, for any x, 

(B IIzS (x)B IIzS (x))^ = _i. 

To continue, we now use symmetries at the level of 
correlation functions to give a more complete treatment 
of the baryon dispersion curves than that appearing in 
Theorem 4 of Ref. 5|. Besides the ordinary time re- 
versal T, charge conjugation C and parity V, which can 
be implemented as unitary (anti-unitary for time rever- 
sal) operators acting on the physical Hilbert space H, we 
also have a time reflection symmetry T as well as the 
combination CV T. 

We summarize these symmetries as follows: 

• Time Reflection T: The ip fields transform as 

1p a (x) -> Aaplpp (~X°, x), $ a (x) -> ij>p(-Xp , X)B 0a 

where A = B = = ^ 2 ) , f(g xy ) - 

f(9sy), with z = (-z°,z); 

• Charge Conjugation C: ip a ( x ) ~ * Af) a , 
i> a {x) -► B a pil> (x), A = -B = B- 1 = 

• Parity V: ip a {x) -> A a pipp(x°, -x), rj) a (x) -> 
ipfs(x°,-x)B0 a where A = B = A -1 = 7 , 
/(Sa^/) -> f(9sy), with z = (z°, -z); 

• Rotation r 3 of tt/2 about e 3 : ipa(x) — > 

A af3 '4> l 3(x ,X 2 ,-X 1 ,X 3 ), $a{x) -> 

tp0(x°,x 2 ,—x 1 ,x 3 )B/3 a where A = _B _1 = 



diag(e-W4, e W4 ;e -W4 )e W4), /( 5x!/ ) _» 
with 2 = (z°, z 2 , — z 1 , z 3 z); 

• Time Reversal T: ip a ( x ) ~ * i^)Ap a , 
^ Q (x) B aP il> fi (-x°,g), A = B = B- 1 = 7 , 

/G?*i/) -> [/(&*)]*• 

The above symmetries are taken to be order preserving, 
except for C and T which are order reversing. For all of 
them, except time reversal, the field average equals the 
transformed field average; for time reversal the trans- 
formed field average is the complex conjugate of the field 
average. 

Applying the VC T symmetry, we find the following 
relations, after making a half-integer shift in the temporal 
coordinate of the lattice, 

(p±(u)p T (v)) = , (p + (u)p + (v)) = (p_(u)p-(v)) , 

(B.5) 

and the same for the n's. As shown in Sec- 
tion we also have (p±(u)n T (v)) = = 
(n±(u)p T (v)) , etc. Setting (f> s (u) = B| 3 s (u), we 
also have (03 (tt)0_| (v)) — (0i(u)0_ 1 (v)) = 0, 
(4> T 3(u)j) T i(v)) = -{(f>±i(u)<i)±z(v)), (0_|(u)0i(u)) = 
(0_i(u)0§(V)), (0|(u)0|(v)) = (<j>_a(u)4>_a(v)), and 
(0a(u)0i (u)) — (4>_i(u)cf>_i(v)) and, using r 3 , for x r = 

( ,-y>0 ™2 ™1 -y>3 \ 

(0|(u)0i(?;)) = -£(0|(u r )0i(?;)) , 

(0|(u)0_i(v)) = -(0|(O0_i(v)) , 

(0§(u)0_§O>)) = ^|K)i|W) , ^ ^ 

(0i (u)4>_i (v)) = -i(<j)i(u r )4>_i(v)) , 
(0i (u)0_|(v)) = -(0i (w r )0_|(w)) , 

(0 s (u)0 s («)) = (<j) s (u r )<j) s (v)) . 

The relations after Eq. (|B.5(I carry over to the two- 
baryon correlation G. 

Using Eq. (|B.5(I . the proton and neutron two-point 
functions are diagonal in spin, at all orders in n. For 
spin 3/2, label the basis by 3/2, -3/2, 1/2, and -1/2, or, 
simply, 1, 2, 3, and 4, respectively. With this labeling, 
and remarking that time reversal T parity V and the 
spectral representation ensures self-adjointness, the two- 
point function matrix has the structure 

/a c d \ 

a d —c I 

c d b I ' 
U -c b J 

with multiplicity two eigenvalues 

M± = {(& + a) ± l(b - a) 2 + 4(|c| 2 + |d| 2 )] x / 2 }/2 . 

The corresponding eigenvectors are (c, d, /x-t — a,0) and 
(m± — b,0,c,d). Note that these two linearly indepen- 
dent vectors have a singular limit when c, d — > 0, where 
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they become identical. [This limit corresponds to taking 
f 13 , f 14 -> in Eq. (|FTH|) below.] 

The same matrix structure is satisfied by its inverse 
matrix and, consequently, it also holds for the Fourier 
transforms G ss '(p° = ix,P) an d F ss >(p° = X S 1. 

Using the r 3 symmetry given above, shows that these 
matrices are diagonal for p = as well as for p = (p 1 = 
0,p 2 = 0,p 3 ), which is used in Ref. [f| to determine 
the one-baryon mass spectrum and to show that no mass 
splitting occurs up to and including order k 6 . For p ^ 0, 
recall that the one-baryon dispersion curves w(p) are the 
solutions of detr(p° = iw(p),p) = 0. From the above 
structure, omitting the p dependence, the determinant 
factorizes as 



detf = [A+A_]^ 



(B.7) 



which gives two by two identical dispersion curves asso- 
ciated to the zeroes of 



2 l 



4 l 



a± = -[r n + r 33 ] ± \l 7 \Tu - r 33 ] 2 + |r 13 | 2 + |r 14 | 2 . 

(B.8) 

For p = 0, A + = Tn and A_ = T 33 , and if they are not 
equal, mass splitting occurs. 

We now show how to find the solution for A + = 0, the 
treatment for other solution being similar. Suppressing 
the k dependence and putting c 3 (p) = c 3 X^/=i 2 cosp 7 , 
for c 3 = —1/8 and c 3 o = — 1 (see Ref. ) , we can write 



r siS2 (p)= [-1 - c 3 (p)n 3 - c 30 K 3 (e^ 



')] 



x5 SlS2 +T,' n > 0tS T SlS2 (n,x)e 

x [tf B0 + (i_ ( y B0 )(e*i' » + e-* ' l )] > 

where ^ means all terms of order k 6 or higher are main- 
tained, and we have used F(x°,x) — T(— x°,x) which 
follows from G(x° , x) = G(—x°, x) by the spectral repre- 
sentation of Eq. 13.1111 . 

To bring the solution in p° to the origin as K \ 0, as 
in Ref. [5j, we make a nonlinear transformation by in- 
troducing the auxiliary function H(w,k) = H(w,K,p) 
and a variable w such that H(w = — 1 — c 3 (p)« 3 — 
c 3 oK 3 e _lp , k) = A+. H(w, k) is then defined by 



H(w, k) = w + 



C30K 



1 ip" 



where 



1 + w + c 3 (p)K 3 ' 2 
x (r n + T 33 )(n, x)e-^- 3 {5 n0 + (1 - S n0 ) 

(B.9) 



means that we only retain terms of order k 6 or higher 
in the sum, and that we have used the short distance 
behavior of r siS2 as given in Theorem 2 of Ref. [| [the 
misprinted values of c 3 and c 3 o in this reference are cor- 
rected above]. 

Before we continue the analysis, we give the following 
short distance behavior that extends the results of The- 
orem 2 of Ref. [5( to a higher order in k (cq here is a 
K-independent computable constants): 



-4c 6 k 6 ,x = (0,1,1,0); 

s, s, = § { 2c 6 k 6 7 x= (0,1,0,1), 

(0,0,1,1); 



si = s 2 



Gs lS2 (x) = <J ( 4c 6 K 6 ,x= (0,1,1,0); 

-2c 6 k 6 ,x = (0,1,0,1), 
(0,0,1,1); 

, Ok 6 , si ^ s 2 • 

(B.10) 

We also find the following k 9 contributions (eg, eg ± being 
nonzero constants): 



G SlS2 (x) 



c 9 k j , x = (0, 1, 1, 1), si = S2 , 
cg,± k 9 , x = (0, 1, 1, 1), si 



±7 



(BT1) 

Note that there are 0(k 8 ) contributions to Ss 1 s 2 ( a; — 
0) arising from two square loops starting and ending at 
x = 0, but with opposite orientations. They give equal 
diagonal contributions and there are no off-diagonal con- 
tributions using the symmetry properties of Eq. I|B.6() . 

If any two of the terms in the square root are zero, 
the square root in Eq. I|B.9|I disappears and, since by 
Eq. I|B.9(1 and the decay properties of T obtained with 
the hyperplane method leads to a joint analyticity in k, w 
and each p component, we can apply the analytic implicit 
function theorem to find an analytic w{p 1 n) such that 
H(w(k,p), K,p) = 0, for small n. In fact, the analytic 
implicit function w(p, k) becomes explicit and admits the 
contour integral representation 



w(j>, «) 



1 

2~Tri 



\w\-- 



w (dH/dw)(w, K,p) 
H(w,K,p) 



dw . 



where the contour can be taken as a small positively ori- 
ented circle (see, for instance, Ref. ^3)- This represen- 
tation is used to obtain an explicit formula for the Taylor 
series coefficients in Ref. |48| . 

Concerning the proton and the neutron, we know that 
their dispersion curves are identical by isospin symme- 
try. Moreover, as the VC T symmetry [see Eq. (|B.5|) ] 
ensures F SlS2 (p) is diagonal (even at p ^ 0) and propor- 
tional of the identity, the proton dispersion curves are 
also identical for their two spin states. The resulting pro- 
ton and neutron dispersion curves are given by, recalling 



Ci = \ E">o.x( r 33 ~ r u )(n, x)e-^- s {5 n0 + (1 - 5 n0 ) that Pt = 2 Ei=i(l " cosp*) as in Eq. El, 



v f 1 \n!( l+w+c 3 (p)K 3 \ n , / C3QK 3 \n]^ 

X I L ) LV > + ll + U, + C 3 (p)K ;1 ) J/' 

and similarly for £ 2 and C3 with the factor [r 33 — TnJ/2 
replaced by Ti 3 and Fi4, respectively. In the above, j^" 



w{p) ee w(p,k) = [-31hk- 3k 3 /4+p 2 k 3 /8] +r 3 (p) , 

(B-12)^ 

where r(p) is of C(k 6 ) and is jointly analytic for each p 
component, and n. As easily seen, w(p) in Eq. (|B . 1 2|> 
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and increases with |p| , as 



is convex for small |p| and 
physically expected. 
As seen in Eq. (E 



TTJ) , we have non-diagonal contribu- 
tions for G SlS2 , |si| ^ | S2 j - The square root in Eq. (IB. 91) 



survives. Using r% symmetry through Eq. ljB.6l) and V, 
we find the following 0(> k 6 ) contributions to G Sl s 2 (p)' 

G| 3 (p) oc (—2 cosp 1 cosp 2 -I- cosp 1 cosp 3 , 

+ cosp 2 cosp 3 ) ft 6 , 

G 1 1 (p) oc (2 cos p 1 cos p 2 — cos p 1 cos p 3 , 

— cosp 2 cosp 3 ) k 6 , (B.13) 

Gi i (p) oc [(1 — i) sinp 1 cosp 2 

— (1 + i) cosp 1 sinp 2 ] sin p 3 K 9 , 

G| _i (p) oc sinp 1 sinp 2 cosp 3 K 9 , 

where we observe that all of these contributions do vanish 
at p = 0. 

We have similar expressions for fy. Hence, !T| |(p) — 

f^(p) 7^ 0, and we cannot successfully apply the same 
method as before, as the k analyticity region shrinks to 
zero as p \ 0. Further analysis is then required. 

Using the above results, without using analyticity, we 
can however obtain the remaining approximate disper- 
sion curves, but with less knowledge about their behav- 
ior. Particularly, we only have a rough bound on r(p). 
To do this job, with h 2 denoting the square root term in 
Eq. (|B.9|) . we can write 



H 



w 



hi + h 2 ■ 



(B.14) 



where, uniformly in p, we have \h±\ < cik 6 and \h 2 \ < 
c 2 k 6 , for some positive constants C12, by the above and 
the bounds of Theorem 2 of Ref. yj . [The h 2 bound is 
due to the k 6 term in Ci , which goes to zero as p \ 0, 
but we can possibly have a term of order 0{> k 10 ) which 
does not go to zero as p \ 0. Such a term leads to a mass 
splitting.] 

Writing H = in Eq. (|B.14|) as 

to = -hi(w, n,p) - h 2 (w, K,p), \hi + h 2 \ < en 6 , 

we plot the left and right sides as functions of w. For 
\w\ < (1 + a)cK 6 , for some a > 0, we have at least one 
solution of H = 0. Similar considerations apply in work- 
ing out the zero determinant solution for A_ in Eq. i|B.7[) . 
From Ref. ,5], a Rouche theorem argument shows that 
there are exactly four solutions of detr = 0, so that there 
is exactly one solution for each A±, which are presented 
in Eq. ETT5) |. 

We now give an explicit construction of the dispersion 
curve in the non-diagonal case. In H{w, n,p), it is conve- 
nient to make another nonlinear variable transformation 
from wtoiiEiii + fi(w, k) and, by the analytic implicit 
function theorem, L(u,w) = w + w + /i(w,k) — u = 
yields the analytic function w = r(u, k) such that 



L(r(u, k),u) = for \w\, \k\ small. The condition for 
the existence of r is satisfied as 

(dL/dw){w, u) = 1 + {df 1 /dw)(w, k) ^ , 

for \w\, \u\, \k\ small, since |/i| < cik 6 and |<9/i/<9u;| < 
c 2 n 6 . Furthermore, (C 2 + IC2| 2 + |C3| 2 )(w = r(u, n), K,f) is 
jointly analytic in u, k, p J and we write it as ^ n a n u n . 

Thus, we can write the dispersion curve determining 
equation H = as the fixed point equation 



E 



a n u 



1/2 



/(«)• 



We obtain the solution u by iteration. The (£ — l)th 
iteration gives 



Sn 2 a ™2 

mi 1/2 



Assuming convergence to the fixed point u, for the ini- 
tial value u = 0, we take the limit i — > 00 in the above, 
with ai replaced by ao Seo, to get the solution u. 

Additionally, by adapting the subtraction method of 
Ref. t 5], it can be shown that the isolated dispersion 
curves (the upper gap property) associated with the two- 
baryon functions of the above 20 particles, with asymp- 
totic mass — 31nK, are the only particles in the full space 
of states generated by an odd number of ip fields, up to 
the meson-baryon threshold w —5 In k (the meson asymp- 
totic masses are ~ — 2\uk). 

To close this Appendix, we point out that for free lat- 
tice Fermi field the matrix structure of the Fourier trans- 
form of the two-point function is the same as in the con- 
tinuum, i.e. p^7'' + ml 2 in the lower subspace. 

From this explicit relation, we see that there is 110 spin 
flip and that the diagonal elements are the same. An- 
other way to see this is to use the VCT symmetry which 
gives (^3^4) = (^4^3) = and (^3^3) = (^4'04>< More- 
over, in the interacting case with a rotationally invariant 
interaction, such as in the U(A^)-invariant Gross-Neveu 
model, the same holds using VCT symmetry. 

However, for the three-quark composite local baryon 
fields, the VCT symmetry, or even the VC T symmetry, 
are not enough to conclude diagonality of the two-point 
function in the spin, except for protons and neutrons. As 
we have seen, it is not diagonal in the spin 3/2 case. 

In the continuum case, relativistic invariance allows 
us to pass to the rest frame, and the two-point function 
transforms according to the spatial rotation group SU(2). 
There is no (spin) mass splitting and there is just one 
dispersion curve, the relativistic one. 

APPENDIX C: Coincident point four- baryon cor- 
relation and linear dependence relation 

Here we obtain the decomposition of Mq in the total 



21 



spin basis, and show the pointwise linear dependence re- 
lation given in Section 4.3 [see Eq. H4.27fl ]. 

At k = 0, using the total spin spatial symmetry and 
orthogonality properties, the following block decomposi- 
tion as given in Eq. I|4.25|) is 

M <°> (0) = Afj® (0) © M[ 0) (0) © M} 0) (0) © M b (0) (0) . 

It is convenient to use the shorthand notation bi zS = 

B(j} IrS . M fe ^(0) is diagonal by the total spin orthogo- 
nality relation and, by the lowering relations of the z- 
component of total spin, all of its elements are equal. 
The first diagonal element is 

,,7-3/2 3/2 ^3/2 3/2 > (01 a iu u 
\ 1 T=3T,=3 1 T=3T z =3l ~ 4 \°3/2 3/20-3/2 3/2 

X &3/23/2^-3/23/2) (0) = ~4 . 

The elements of the 2 x 2 matrix M 1 (0) (0), with 5=1, 
S z = l,T = l,T z = 1, are given by 

/in = 2(p + n + p + n+) (0) , ^22 = (xx) {0) , 
M21 = M12 = (p+n + XX) {0) , 



/2 1/2 



V3/1 Oi /23/ 2 
1/21/20-1/21/2 - \/3/10 6i/2_i/2 



where x = y /3/10 6 3 / 23 /2 6_ 3 / 2 _i/2 

0-3/21/2 + \/3/10 3 / 2 _i/2 &_ 3 / 23 / 2 ' 

0-1/2-1/2 + V 2 / 
0-1/23/2- 

For the sake of checking the calculations, we compute 
the term of the four-baryon correlation appearing in fi^j 
with 6 fields which have the prefactor 1 since their values 
are integers. We denote them by B u and call them un- 
normalized fields, and label the unnormalized products 

n U U UU UU y-a IK Ltl LU 1,11 1M 

P+''_l_, 3 3 U -3 -1] l<3 1 Ml] "2^1 3 1 "1 3 O-l -1 ; 

22 22 22 22 22 22 22 22 

b\ i i , 6" _i 6^1 3 appearing above by 1,2, . . ., 7, 

22 2_ 22222 

respectively. Here, the corresponding n = averages are 
denoted by (i,j), where the second spot corresponds to 
the product of two barred fields. By isospin flip sym- 
metry, we have the following equalities: (1,2) = (1,4), 
(1,5) = (1,7), (2,2) = (4,4), (2,6) = (4,6), (3,5) = 
(3, 7), and (5, 5) = (5, 7). Thus, all we need to calculate 
are the elements listed below with their respective values: 



(1,1) 


= -360 , 


(1,2) 


= 144, 


(1,3) = 


-96 


(1,5) 


= -144, 


(2,6) 


= 288, 


(3,3) = 


-144 


(3,5) 


= 96, 


(3,6) 


= 240, 


(4,7) = 


144 


(5,5) 


= -240 , 


(5,7) 


= -192 , 


(6,6) = 


-1296 



Reinstating the normalization prefactors and substi- 
tuting in fiij results in the matrix M^°\o) of Eq. H4.26fl . 

We now turn to the linear dependence relation of Eq. 
KTfy . To prove that 



V5 



S=1,S* 



T=l,T, 



o, 



for S z 
for S z 



T, 



— 1,0,1, it is enough to show it holds 
1 and apply the spin lowering operator 



on S z and T z . For S z = T z = 1, using the definition 
of the T fields in Section 4.1 [see Eq. (|4.2U|) ]. we write 

Ts=i i + =r^r=i l as the sum, with from now on 6r s = 



B(I=3/2)J.a, 5 = ±3/2, ±1/2, 



03 -1 6-3 3 

2 2 2 2 



6336-3-I " 

2 2 2 2 

6136-I -1 
2 2 2 2 

bi -1 b-i 3 

22 22 



1 63 16-3 1 



+ \ h 6116-1 1 

\ z 2 2 ~2~ 2 



A typical term in the above is labeled by a product of 
six -0's, with indices from one of the following 12 sets of 
indices: (1) 1 - -; (2) 1 - +; (3) 1 + -; (4) 1 + +; (5) 
2 - -; (6) 2 - +; (7) 2 + -; (8) 2 + +; (9) 3 - -; (10) 
3-+; (11) 3 + -; (12) 3 + +. 

There can be no repeats of any of the twelve indices in 
the six factors, by Pauli exclusion. We are left with 21 
linearly independent products of six ^'s with no repeated 
indices, and the coefficients of each of these linearly in- 
dependent elements sum to zero. 

APPENDIX D: Decay of the Bethe-Salpeter kernel 

Using the hyperplane decoupling method, as described in 
Appendix B (see also Ref. [5[ for more details), we obtain 
improved temporal decay for the kernel K, as compared 
to that of M and AIq. For now, we use the original and 
unmodified M, and we will treat the modified M = M' 
at the end of the Appendix. 

The temporal decay follows from a product structure 
of the sixth derivatives of the four-point functions M and 
Mo, which in turn follows from the product structure of 
the sixth k p derivative of the basic four-point functions 
D and Dq defined in Section 4.2. Recall that D, in the 
individual spin basis is defined in Eq. 14.22J1 and Dq is 
obtained from D by erroneously applying Wick's theo- 
rem to the baryon fields in D. In a similar way, Mq 
is obtained from M. In what follows, we suppress the 
collective index subscripts from D and Do. The prod- 
uct structure extends from D (Dq) to M (Mq) as it is 
preserved under real orthogonal transformations. 

We show that, for = < p < s 3 = and with 
w, W\, W2 € Z 3 , 

D^( Xl ,X 2 ,X 3 ,X 4 ,) = -lEtfi^a 

x D ( -°\x 1 ,x 2 , (p, wi), (p, w 2 )) 

x D(°\(p + l,w 1 ),( P + l,w 2 ),x 3 ,x 4 ) (Dl) 

-l£tf D( \xx,x 2 ,(p,w),(p,w)) 

x £>(°)((p + 1, w), (p + 1, w), x 3 ,x 4 ) , 



and 



D ( 6) (x 1 ,X2,X 3l X 4 ,) = -^T,w u w 2 

x D^ ) (x 1 ,x 2 , (p,wi), (p,w 2 )) 

x D Q > \(P+ 1, ^i), (p+ l,w 2 ),x 3 ,x 4 ) 



(D2) 
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and similarly for x\ > Xg. We write these results 
schematically as 



D W = -l D W ® D (o) _ l D (o) @ D (0) 
2 4 



(D3) 



and 



(0) 



)(0) 



(D4) 



©" and "o" 



using the convolution-like operations 
given by Eqs. JHB and (|D2)l . 

We now show how this structure gives the improved 
temporal decay for K = Dq 1 — D~ x . Take x\ < p < 
X3, for example, and let A (A ) denote the inverse of D 
(Do). It is easily seen that K^> = 0, n = 0, 1, . . . , 5, 
which follows from the fact that the same holds for D 
and Do- This is seen by expanding the numerators and 
using imbalance of fermions. We have, using Eqs. I)D3|1 
and HD4|I . 



K<& = A^D^A^ -A(°)D( 6 )A(°) 



= +Ar(-iDr°4 o) )A o) 

_A(o)(_i£)(o) ®d(o))a(o) 
-A(°)(-i J D( )®D(°))A( ) -0, 



(D5) 



for x^+l < X3. As this holds for each hyperplane between 
and x°, the above leads to a k 6 k 7 ^ -2 ! -1 ) decay, for 
x§ > x? + 1 and a decay k 6 , for X3 = x\ + 1. The ft 7 can 
be further improved to k 8 taking into account fermion 
imbalance (or gauge integration) in the calculation of the 
seventh k p derivative. 

We return to show how Eqs. (|Dlfl and (|D2() are ob- 
tained. For x\ = x§ < x§ = x°, 



D, 



^= -(B tl ( Xl 



3 (x 3 ))( 3 ) (B e2 (x 2 )B ti (x 4 ))(V 
B ei ( Xl )B ei (x 4 ))^ (B t2 (x 2 )B e3 (x 3 ))W , 



where the £'s are collective indices. 

Using Eq. (|B.3|) . and regrouping the terms, Eq. (|D2|> 
follows. 

To show Eq. (|D1(I . in addition to the gauge integral X3 
for three coincident bonds of Eq. 1|B.1(I . we need to com- 
pute the gauge integral T§ of six coincident bonds (same 
orientation!). To compute X 6 , we use the generating func- 
tion given in Refs. and differentiate with respect to 
the external sources. It is given by, with U (g) = g and 



U-\g)=g-\ 



(D6) 



50161 9a 2 b 2 g as b 3 9 A. 9 A. 9 A. dn(g) 



3!4! ^ 6ai a2a3ebl b2 63 £a 4 a 5 a 6 e b 4 b 5 b 6 - 
£01 o 2 o 4 £&i 6 2 64 £a 3 a 5 o 6 £63 65 6 6 + 
£aia 2 05 £616265 £a 3 o 4 o 6 e b3 6466 
e aid2a 6 ^bib 2 b e £-0,3(14(15 £63 63 65 + 
£010304 £6163 63 £02 0506 £62 65 66 
£010305 ^616365 £020406^626466 4" 
£ai a 3 o 6 £61 63 6 6 £a 2 04 05 £6 2 64 &5 

£010405 £616465 ^020306 ^62 63 ^6 ~t~ 

£010406 £616466 £020305 £626365 ~t~ 
£oio 5 o 6 £6i656 6 £a 2 0304 £626364] • 



By performing appropriate (index) contractions inXg, we 
verify that X4, X3 and X2 are obtained, where (see Refs. 

USES) 



1 



50161 9 Abo ^(9) = o <W 2 ^261 j 



(D7) 



X3 is given in Eq. IjB.lf) and 

X 4 = /.9o 1 6 1 5q2 1 6 2 .9o 363 3o4 1 6 4 ^(5) 

= g [<5ai6 2 < 5a364^6io 2 ^6304 + ( a 2 ^ a 4, ^2 ^ &4)J 

— 24 [^0162^0364^6104 ^6 3 a 2 + ( a 2 ^ 04, 6 2 ^ 64)] 



(D8) 



Eq. (|D6|) can be written in a compact notation as 



Ie = r3!4! ^ S * £ s> 



where T 3 is the set of the ten disjoint partitions of T = 
{1,2,3,4,5,6} = iJU 5, = 3 = and for R = 
fc}, B fl = e £6i6,6 fc , where the elements in _Er 
and fig are written in increasing order. 

£>( 6 ) is obtained by writing D — N/Z and expand- 
ing the exponential in the numerator of N . In this way, 



r(6) 



V 1 

1 *n 



(6) 



£>(«) = A^( 6 )/Z(°), where = M 6) + jvj' 

arises from the terms where there are two nonoverlapping 
bonds, indexed by w\ and W2, each bond with three pos- 
itively oriented gauge fields. After performing the gauge 
integrals using the formula for X3 of Eq. (|B.1() . we get 

»t(6) 1 yv 

n ~ 2xgT £010203 £o 4 0506£6i6263 £6 4 6 5 6 6 2-*,w 1 ^w 2 ,{a i },{fi} 

X (B ll (x 1 )B t2 (x2){^a 1 a 1 f 1 ^Pa 2 0:2/2 ^0303/3 ) (P. ^l) 

X (V>o 4 a 4 /4 ^050-5/5 ^a 6 a 6 / 6 )(P,^2)) (0) 

x ((VW/i i>b 2a2 f subsets f 3 )(p+ 

)(p + 1, ilia) 

x4(i3)4(^)) (0) , 
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where the numerical factor comes from ( — l/2) 6 x(l/6!)x 
20 x (-2) 6 x (1/2) x (1/6) 2 ; 20 is a combinatorial factor, 
(— 2) 6 comes from the six T e 's and (1/6) 2 coming from 
the two integrals J3. We observe that the sums over the 
Qii's are restricted to lower components a = 3,4 due to 
the structure of T e . 

Resumming over the intermediate states, after tak- 
ing into account multiplicities and normalizations of the 
baryon fields, gives the first term of Eq. (|D1|I . 

The contribution to N%, for which there are overlap- 
ping bonds, after performing the 2q gauge integral, is 

~~ 4x6 4 x 10 ^w,{ai},{fi} [ e a 1 a 2 a 3 £b 1 b2b 3 £a4,a s a 6 £b 4 ,b 5 b 6 

+ nine remaining terms] (B^ {x\)Bi 2 (x2) 

X (^0101/1 ■ ■ ){p,w)r 0) 

x ((iphah ■ ■ ■'0h 6 a 6 / 6 )(p+ l,w)Be 3 (x 3 )B l4 (x4)) iQ) , 

where the numerical factor comes from ( — 1 / 2 ) 6 x ( 1 / 6 ! ) x 
(— 2) 6 x [2/(3! x 4!)]; the term in square brackets is the 
contribution of the integral. Again, resumming over 
the intermediate states gives the second term of Eq. 
(|D1(I . In both cases, a factor 6~ 4 is absorbed in the 
normalization process. 

By considering hyperplanes in each spatial direc- 
tion, and additional complex parameters associated with 
them, spatial decay of the B-S kernel can also be obtained 
as in Ref. 0. 

We now describe how the product structure of Eqs. 
(|D3)l and (fD4)l extend to and M (6) . First, we take 

isospin C-G linear combinations of the external fields 
present in D^ & \ expressed in the individual spin and 
isospin basis, to get total isospin zero states. Then, again 
using the inverse C-G coefficients, we express the inter- 
mediate states in terms of states of fixed total isospin 
I and its z-components I z . At this point, the isospin 
orthogonality properties are used to get only zero total 
isospin intermediate states. In this way, we arrive at the 
following product structure formula for 

M< 6 > = -Im<® © M<°> - - A M^ ® M<°> , 

where we remember that we are still in the individual 
spin basis. We can pass to the total isospin zero, total 
spin basis by taking spin C-G linear combinations. The 
product structure is maintained as it is invariant under 
real orthogonal transformations. A similar argument ap- 
plies to 

We are now ready to explain and motivate the adoption 
of a /i-modified M = M' as introduced in Section 4.3. For 
this, we consider the following leading contributions to 
the unmodified K(x\, X2, X3, X4) and the corresponding 
contributions to the generalized potential K(£,rj,k°): i) 
the k° order, K^ ' at coincident points, and correspond- 
ing space-range zero, energy independent constant poten- 
tial K(0, 0, fc°); ii) the k 6 contribution with time-range 
one and space-range zero K^ 6 \0, 0, e°, e°), with energy- 
dependent zero space-range potential K(0, 0, k° = i(2fh— 



e)) of order K 2 e _c ; and, in) the k 2 quasimeson exchange 
n 2 contribution of space-range one K^ 2 \0, 0, e J , e 7 ), j — 
1,2,3, with space-range one energy-independent poten- 
tial K(0, e j , fc°) of order k 2 . 

So, the resulting potential is seen to be energy- 
dependent. We make our modification to M and then 
K so that only an energy-independent potential arises 
which is more in keeping with our usual intuition for un- 
derstanding the occurrence of bound states. 

As the fc°-dependent part of the effective potential is 
associated with a temporal distance-one k 6 contribution 
to K , wc modify M so that this contribution is smaller. 
This happens because M is modified so that M^ 6 ' has the 
same product structure as Mq . By doing so, the decay 
of the modified K is improved for temporal distance one. 
In contrast to the unmodified M, for our modified M = 
M', we get a k s ^ x ^~ x ^ decay for any distance \x® — x\\ > 
1. This comes about as the modified also has a 

product structure so that the modified K (x\, X2, £3, X4) 
in Eq. (|4.28l) has a zero sixth derivative = 0, for 

x\ = X2 < £3 = since the last two terms on the r.h.s. 
of Eq. (|D5(1 have now the same coefficient as the first 
term. 

There is still a possibility of having a k 8 contribution 
to the modified K (x%, X2, £3, X4) of temporal distance one 
leading to an energy-dependent potential of order k 2 , 
which would be of the same order of the quasi-meson 
exchange potential. These k 8 contributions arise from 
M and Mq with temporal distance one, but vanish due 
to imbalance of fcrmion components. 

With this modification, the space-range zero potential 
is zero and the dominant contribution is the quasimeson 
space-range one exchange potential. 

APPENDIX E: Determination of the space range- 
one exchange potential 

In this Appendix, we determine the quasi-meson ex- 
change potential 20 x 20 matrix of Section 4.3 in the 

individual spin basis, i.e. ^ 2) S2S3S4 , v^ aataU , v£l 2SaSA , 
(2) 

and v\ t 2 uti- We give reduced formulas which we use 
for computation, and present the derivation explicitly for 
vi 2 \ 2 s 3 s4 given in Eq. (|4.31|) . the elements of which, in 
the total spin basis, are given below Eq. I|4.32|l . The 
other elements are obtained in an analogous manner. 

In Eq. H4.31fl . consider the K — averages expressed 
as a sum of all contractions using Wick's theorem. If the 
explicit tjjij) contract in the first average then cti = [3% 
and r^ i/3i r~ 2 e Qi = 0; the same for the second average. 

Thus, only contractions involving one explicit internal ip 
and one quark field of B's contribute. Also, recalling 
that the ^'s on the B's only have lower spin indices, this 
forces the spin index of the explicit ip's to be also lower in 
order not to get a zero contribution. Next, noting that in 
the subspace of lower components the spin matrix r ±e 
reduces to minus the identity, and after using the k = 
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isospin flip symmetry in the second factor, we can write 



(2) _ 3 3 3 3 

V Sl ,s 2 ,s 3: s 4 — C I 1 -I 1 C L 3 -I 3 X 

Il,h 



(Bi Ilsl Bi l3S3 : i> aafl ^ a pf 2 :) (0) X 
(Bi IlS2 Bi hSi : ij b p_ h ij) ba _ h :>(°) 



(El) 



with all fields at a same site and where the Wick or- 
der : : forbids contraction between the explicit internal 
ip's. To implement the exclusion of this contraction, we 
note that up to a global sign and omitting the detailed 
field indices and summations we have (BBipip)^ ! ~ 

5, ...,8, such that the Wick 



det[5j. 



1, ...,4 and j 



order in : V^V's : is equivalent to replacing Sys by zero in 
the previous fermionic Gram determinant detftfy]. 

To continue, we use symmetries and isospin and spin 
restrictions to reduce the number of u' 2 )'s to be cal- 
culated. Fixing a matrix element in v^ 2 \ we consider 
the terms in the sums labeled by {Iih), (I1J3), (Jih) 
or (J1J3). Isospin flip symmetry gives equality of the 
(J1J3) and {—J 1 — J3) terms, etc., where we take into 
account the properties of the C-G coefficients. Isospin 
conservation forces the restrictions I\ + f% = I3 + /1 => 
\h - h\ < 1, h + h = Jz + h => \h - J 3 \ < 1, 
J\ + f-2 = J3 + fx => \Ji — J 3 \ < 1. In addition, we 
have spin sum restrictions si — S3 = S4 — s% = a — /3, 
\si - s 3 1 < 1; si - t 3 = ti - s 2 = a - P, |si - t 3 \ < 1; 
ti — t?, — <4 — i 2 = a — (3, |ii — £3] < 1. Also, for k = 0, 
conjugation (denoted by Co) is a symmetry. Using Cq 
gives the equality of the (/1/3) and (/3/1) ((Ji J 3) and 
{J3J1)) product of summed expectations if s± = S3 and 
•?2 = S4 (*i = ^3 and £2 = £4)) such as diagonal elements. 
In this way, we show that many of the terms are zero. 
We can also, after applying Cq to both factors, change 
/1 <-> /2 <-> -/a to get equality of the (7i/ 3 ) and 

(J3J1) terms if si = S4 and S2 = 53; similarly for the J's. 
We can also impose the prohibited contraction condition 
to show that certain elements are zero. 

In addition to the fact that the 20 x 20 potential matrix 
is symmetric there are other symmetries. By making the 
change of labeling a <->• @, fx <-> ~f 2 , f-2 <-> the 
following symmetry properties are obtained 



.,(2) 



- J 2 ) 



V 



(2) 

; 

si s 
(2) 



7' 



(2) 



(2) 

t±t2t 3 t4 u t2t\t&t 3 



Also, the spin flip change a <-» —a, /3 <-> — /3 gives the 
equality of the v( 2 ^'s for Sk — > — s^ — > — for all j's 
and fc's. 

Thus, there is symmetry in the upper left 4x4 (lower 
right 16 x 16) matrix about the secondary diagonal of the 
submatrix. 

Using each symmetry alone or coupling them explains 
the repeated values occurring in the matrix with the 
following two exceptions: the equality of the (6 9) and 
(811) elements and the equality between (1111) and 



(11 14). By spin flip at k = 0, we see that (1 1) = (4 4), 
(2 2) - (3 3), (17) = (418), (110) = (415), (113) - 
(412), (2 8) = (317), (211) - (314), (214) = (311); 
by spin flip at k = and 1 <-> 2, 3 <~+ 4 exchanges 
(113) = (418), (7 7) = (1212); by 1 <-► 2, 3 <-► 4 ex- 
changes (6 6) = (9 9). By applying Cq to both average 
factors, we obtain (710) = (1013). Thus, we only have 
to calculate e.g. the elements (1 1), (22), (23), (17), 
(1 10), (1 13), (28), (211), (214), (55), (66), (69), (77), 
(710), (88), (811), (1010), (11 11) and (11 14). 

The values of the above elements are given below Eq. 
(|4.32l) . As a check in calculating the elements of Eq. 
(|4.31l) . etc, it is useful to give the decomposition of the 
total value before performing the sums in the equivalent 
expression of Eq. I|E1|) . We now give the nonzero values 
for the partial terms (ii/3) or (Ii J3) or (J1J3) for each 
matrix element; the values for (— 1\ —I 3 ) or (— 1\ — J3) or 
( — J \ — J3) are the same. The sum of the values for each 
relevant partial matrix element given below multiplied 
by a factor 2 is the value of the matrix element. The 
values, given here with precision of two digits, are as 
follows. For the matrix element (11), which corresponds 
to si = S2 = 1/2, S3 = S4 = 1/2, in Eq. (|E1|) . we need the 
contributions (/1/3) = (5 |), which is 0.78, and (ii/3) = 
(i — i), which is —0.94. Denoting this contribution as 
(11): {- 2 2 ; 

-0.44, (I -1) -1.39; (17): (If) 



1 )) 0.78, (| - |) -0.94, and using this notation 
for the other elements we find (2 2): (L, . — , - 
i) 0.44; (2 3): (± \) 

-0.41, (ii) -0.27, (i - \) -0.14; (110): (± §) 0.47, 
(i i) 0.31, (i -i) 0.16; (1 13): (± f ) -0.41, (± \) -0.27, 



$ -I) -0.14; (2! 



(i |) 0.71, (ii) 0.47, (|-i) 0.24; 



-0.16; (214): 



^2 2) vj.^-l, y 2 2 > 

ls 0.08; (55): (| ±) -0.75, 

'3 1 j 

(i 3) „ . 50j ( i i) . 67j ( i _i) „ . 67 . (69 ). (| i) _ . 25) 

(i |) -0.25, (i |) 0.33, (i -|) -0.33; (77): (§ i) -0.25, 
(||) -0.25, (ii) 0.33, (i - i) -0.33; (710): (| ±) 



(211): (if) -0.47, (|i) -0.31, (i 
(II) 0.24, (|i) 0.16, (i 2, 

(i |) -0.75, (i i) 1.00, (i -i) -1.00; (66): (| |) -0.50, 
\) -0.50, ( 
2 2") "0.25, ( 2 _> 



(if) "0 2 29, (|i) "0.38," (i 



^2 2> 

-0.29, 

all zero; (811): (§ \) -0.25, (||) 



i) -0.38; 



2 

-0.25 



i) -0.33; (1010): (| \) 



(ii) 0.33, 



-0.42, (if) "0.42, (|i) 
-0.33, 



(i - i) 

V2 2^ 

0.56, (| - i) -0.56; (1111): -0.33, (if) 

(ii) 0.44, (i - |) -0.44; (1114): (| I) -0.33, (if) 
-0.33, (|i) 0.44, (i - i) -0.44. 

The values given below Eq. I|4.32|l are the exact values 
for ii ( 2 ) obtaining by summing up the exact values for 
the above partial decomposition values and multiplying 
by two. 

APPENDIX F: Uncorrelated four-baryon function 
and spectral representation 

Here we show how to obtain a spectral representa- 
tion for M( it i 1 i 2 i 3 i 4 (x±, X2, X3, X4), and its total spin ver- 
sion; the components of which are obtained by sum- 
ming over the appropriate C-G combinations of the 
Fo,iilal 3 i4,{ x U x 2i x 3) x ^)- Our starting point is Eq. 
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-^0,1234 — — G13G24 + G14G 



23 



where we use a shorthand notation for the collective in- 
dices (£, x). 

In lattice relative coordinates and in the equal time 
representation 



Fo,vsu(Z,ri,T) = -G 13 (t + £)G 24 (r + ff) 



(Fl) 



and setting k — in the Fourier transform of r only, 
with dual variable k, our goal is to obtain a spectral 
representation for 



^0,1234 (I*) ff,k°) = / F , 1234 (i,f},r)e 



dr. 



For a; 7^ 0, and using the spectral representation for 
Gij (x) of Eq. (|3.11|) , we obtain 

Gij(x) = - [ f (A )I^Vi e ^ dAoa ^.(A°)dp, 
J-x Jt 3 

(F2) 

where, by parity, a$,ij = ot-p,ij. 

Separating the r° = contribution and substituting 
the spectral representation for the G's into Eq. i|Fl|l . 
after performing the integral over r, we get 

£o,i234(£ V,k°) = (2ir)- 3 f T3 [-G 13 (p)G 24 (p)E + (p) 
+G 14 (p)G 23 (p) S_(p)] dp 

x [-^Aoap,i3(^°) dxmaff t 24(A'°) 

+d A oap.. 14 (A°)d A /oaj ?i 23(A' ) H_(p)] dp, 

(F3) 



where we recall that G(p) = e lp x G(x° = 0, x), and 
where we have 



S±(p) = S±(p, £, 77) = cosp.£ cosp.ff ± sinp.£ sinp.77 . 

To proceed in our deduction of Mo, we make the ap- 
proximation that G is diagonal in spin, and independent 
of the isospin / and spin (which holds up to and including 
0(k 5 ), see Ref. 0), and we denote it simply by G(p), 
with no index; and the same for the measure d\oap. 

Under this hypothesis, carrying out the sum over the 
isospin C-G coefficients we obtain, for the 11 element of 
Mo in the individual spin basis 

M 4 MlS2S3S4 {£,rj,k Q ) = 

(2tt)- 3 J t3 [G(p)} 2 [-6 s1S3 6 S2S4 Z + (p) 
-6 SlS4 6 S2 s 3 E-(p)]dp 
+ (2 7 r) 3 ^ 1 ^ 1 / T3 /(fc ,A«A'0)x 
[-5 SlS3 6 S2Si 3+(jp) - S SlS J S2S3 E-(p)] 
x d\aafi(\°) d\toap(\ )dp . 

( F4 ) 

The 22 element of Mo is obtained in the same way, by 
replacing Sj by tj in the above. The off-diagonal elements 
are zero. 

Finally, we obtain the total spin basis representation 
for Mo by carrying out the sums over the spin C-G coeffi- 
cients and using the explicit values and symmetry prop- 
erties of the C-G coefficients, we obtain the expression 
given in Eq. (|4.36(l . 
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